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Abstract 



We show how to construct non- spherically- symmetric extended 
bodies of uniform density behaving exactly as pointlike masses. These 
"gravitational monopoles" have the following equivalent properties: 

(1) they generate, outside them, a spherically-symmetric gravitational 
potential M/ \ x — xo| ; (h) their interaction energy with an external 
gravitational potential U(x) is -Mf/(xo); and (hi) all their mul- 
tipole moments (of order i > 1) with respect to their center of mass 
O vanish identically. The method applies for any number of space 
dimensions. The free parameters entering the construction are: (1) 
an arbitrary surface S bounding a connected open subset of K 3 ; 

(2) the arbitrary choice of the center of mass O within £1 ; and (3) the 
total volume of the body. An extension of the method allows one to 
construct homogeneous bodies which are gravitationally equivalent (in 
the sense of having exactly the same multipole moments) to any given 
body. 

Though our method generally assumes that the domain Q is boun- 
ded (which leads to bounded monopoles with closed internal cavi- 
ties), it can also generate unbounded monopoles, with exponentially- 
decreasing thicknesses at infinity, and cylindriclike internal cavities 
reachable from infinity. This may be useful for optimizing the shape of 
test masses in high-precision Equivalence Principle experiments, such 
as the planned Satellite Test of the Equivalence Principle (STEP). 
By suppressing the couplings to gravity gradients, one can design, 
with great flexibility in the choice of shapes, differential accelerometers 
made of nested bodies, which are (exactly or exponentially) insensi- 
tive to all external gravitational disturbances. Alternatively, one can 
construct nested bodies of arbitrary densities having identical (or pro- 
portional) sequences of multipole moments, thereby also suppressing 
any differential acceleration caused by external gravity gradients. 

1 INTRODUCTION. 

The Equivalence Principle states that test bodies of different compositions 
must "fall" exactly in the same way in an external gravitational field. This 
fundamental principle, which lies at the basis of Einstein's theory of grav- 
itation, deserves to be tested with the utmost precision available. Testing 
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the Equivalence Principle is also the most sensitive way to search for hypo- 
thetical non-electromagnetic long range forces. Such forces could be due to 
the exchanges of new spin-1 or spin-0 particles (as suggested in some super- 
symmetric or superstring theories 0,0), whose existence might be required 
for a consistent theory of quantum gravity. Because of their composition- 
dependenceQ, these forces would, by superposing their effects to those of 
gravitation, lead to (apparent) violations of the Equivalence Principle. 

Improving upon the results of the classic Eotvos and Dicke experiments, 
recent laboratory experiments have checked the validity of the Equivalence 
Principle at the ~ 3 x 10~ 12 level In addition, the Lunar Laser Ranging 
experiment has established that the Earth and the Moon "fall" towards the 
Sun with equal accelerations, to within rs 5 x 10 -13 |4]]. A proposed satellite 
experiment known as STEP (Satellite Test of the Equivalence Principle) [0], 
presently studied by ESA, CNES and NASA 0, [/J ||, |], aims at improving 
the sensitivity down to the w 10 -17 level. 

The STEP experiment consists essentially of measuring the differential 
accelerations between several pairs of nested cylindrically-symmetric coax- 
ial test masses of different compositions freely falling within a "drag-free" 
satellite in low orbit around the Earth. At the impressive sensitivity level 
of 10~ 17 , one has to worry about many possible disturbances. In particular, 
the coupling to external gravity gradients of the two differently-shaped test 
masses making up one differential accelerometer causes a parasitic differen- 
tial acceleration even if the Equivalence Principle holds exactly. The origin 
of this disturbance is that an extended body is generally not equivalent to 
a pointlike mass in its interaction with an external gravitational field. In- 
stead, the gravitational structure of an extended body can be described by 
the infinite set of its multipole moments, starting with the monopole moment 
- identical to the total mass. Choosing as origin the center of mass of the 
body ensures the vanishing of its dipole moment (i = 1) and allows one to 
describe its extension-dependent interactions with an external gravitational 

4 New interactions due to spin-1 exchanges are generally expected to act on a linear 
combination of baryon and lepton numbers, or more specifically on their difference B — L , 
in grand- unified theories |TJ. The spin-0 exchange forces suggested by string theory are 
expected to act effectively on a linear combination of B , L and of the nuclear electrostatic 
energy Q. New spin-1 or spin-0 induced forces could also act, in addition, on particle spins. 

5 The common vanishing of the dipole moments of two nested masses - i.e. the fact that 
their two centers of mass can be made to coincide - is realized to a high precision in the 
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field in terms of its successive multipole moments Qi of order t > 2 . The 
existence of a moment of order I leads, in the presence of a disturbing mass 
at distance 71, to an extension-dependent acceleration oc Qz/IZ^ 2 . 

These parasitic accelerations are a serious concern in the STEP experi- 
ment. The way this problem is currently addressed is multi- pronged. First, 
while the inner mass M in is shaped as a straight cylinder with vanishing 
quadrupole moment = 0, but non-zero (even) higher moments Q™, 
Qg n , . . . , the outer mass M out is shaped as a "belted cylinder" [0, ||, |TTJ 
having a vanishing quadrupole moment, Q% nt = 0, and a non-zero hexade- 
capole moment chosen so that Q£ ut /M ont = Q^/M m . The latter matching 
condition for the reduced hexadecapole moments Q^/M ensures the can- 
celling of the differential acceleration oc 1/7Z 6 induced by the two Q^s. 
The remaining free parameters are optimized so as to minimize the effects of 
the non- matched higher multipole moments fllT] , [T2[. Second, from the cal- 
culation of the residual sensitivity of differential accelerometers to external 
masses, one deduces constraints on the allowed motion of nearby masses. In 
particular, one has still severe constraints on the amplitude of any "helium 
tide" in the spacecraft's dewar surrounding the test masses. One must then 
design the dewar so as to meet these constraints. 

In this paper, we propose a radical solution to the problem of the sen- 
sitivity to external gravity gradients. We show how to construct extended 
homogeneous bodies of very general shapes^] behaving exactly as pointlike 
masses in any external gravitational field. The existence of such aspherical 



gravitational monopoles has been made plausible by Barrett [ 13 1 wno found 
approximate numerical solutions to this problem. Here we shall prove the 
existence of such objects by providing an explicit method for constructing 
themQ. Our construction allows us to have a clear control of the large flexi- 

STEP experiment by measuring the effect of the Earth's gravity gradient and positioning 
the masses so as to eliminate this effect. 

6 This very large flexibility in the choice of the shapes of our monopoles could make 
them of practical utility in Equivalence Principle tests. The obvious solution of spherically- 
symmetric test masses was considered by the STEP team |^] but rejected because of 
the unsurmountable practical difficulties of building nested spherical test masses with 
independent suspension and sensing systems. 

7 After submitting this work for publication we learned of an unpublished work by 



Barrett 14 , in which he proposed the same construction method (for monopoles with 



boundaries having the topology of a sphere) without, however, developing it and studying 
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bility in the shapes of such homogeneous^ monopoles. 

An extension of our method allows one to construct gravitationally similar 
bodies (of equal or different densities), i.e. bodies having the same infinite 
sequences of reduced multipole moments Qi/M . Two such bodies experience 
exactly the same acceleration in the presence of arbitrary external gravity 
gradients, as soon as their centers of mass are made to coincide. 

As far as we know, the specific problem of the construction of aspherical 
monopoles has not been dealt with in the large literature on potential the- 
ory. (For an entry into this literature see e.g. Ref. fllBfl .) In fact, there are 
theorems in inverse potential theory proving, under certain assumptions, the 
inexistence of such monopoles. In particular, a famous theorem of Novikov 
[ |T7| ] states that, if we know the value p of the supposedly uniform den- 
sity, the exterior potential generated by a starlik^\ open set Q (with known 
center O) determines uniquely that domain. A corollary of this theorem 
is that a homogeneous starlike monopole must be a ball. We shall in fact 
prove below that this inexistence of aspherical homogeneous monopoles holds 
under the much weaker assumption that the open set Q be connected and 
have a connected boundary S = dQ. I.e., aspherical solid homogeneous 
monopoles cannot exist (this excludes, for instance, monopoles in the form 
of solid tori)0. 

This no-go theorem, however, does not exclude aspherical monopoles - 
either simply connected^ or not - having a disconnected boundary £ = 
E ext U Si nt . Indeed, we shall construct such homogeneous aspherical mono- 
poles, which have a closed (in the sense of being unreachable from infinity) 

its consequences. 

8 By contrast, the existence of inhomogeneous aspherical monopoles is not surprising as 
one can smoothly deform, within some sphere, the potential generated by a homogeneous 
ball, and then define the density "generating" this deformed potential (which remains 
proportional to 1/r outside this sphere) by Poisson's equation. 

9 An open set in M. N is said to be starlike with respect to a center O if O contains 
the straight segment joining O to any point of f2 . 

10 After submitting this work for publication we learned from a referee of an unpublished 
work by Barrett jl5| , in which he established this result by a different, more involved 
method. 

n In Ref. |l7| one finds, without proof (on the ground that it is "easy to prove"!), the 
incorrect proposition that a homogeneous simply connected monopole (with C 2 boundary) 
must be a ball. 
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internal cavity. 

Our construction allows us to choose arbitrarily, either the outer bound- 
ary, or the inner boundary (surface of the internal cavity) of the monopole. 
The topology of these boundaries does not need to be that of the sphere S 2 . 
For instance, the outer boundary can be toroidal. In the simplest^ bounded 
monopoles we construct, the center of mass - which is the center of symmetry 
of the generated external gravitational field - can be systematically located 
within the internal cavity. This allows for a second monopole to be placed 
inside the first one, with the same center of mass. But this inner monopole 
is then unreachable from the outside by a continuous path. It means that, 
in practical applications such as STEP, where one needs to control and sense 
the position of nested monopoles, one can only use approximate monopoles. 

Interestingly, we can get practically useful shapes by considering a limit 
where the boundaries are not compact, but extend to infinity. In particular, 
we shall discuss in detail the case (directly relevant for STEP) where both the 
internal and external boundaries asymptotically approach a straight cylinder 
at infinity. The thickness of the corresponding monopole, around this limiting 
straight cylinder, then tends exponentially to zero at infinity. By cutting 
off such a monopole, one can construct a finite, cylinderlike object whose 
couplings with gravity gradients are suppressed with an exponentially good 
precision. 

After formulating the problem in technical terms and proving the inexis- 
tence of "solid" aspherical homogeneous monopoles in Section 2, we present 
our method for constructing "hollow" aspherical monopoles in Section 3. We 
also indicate how, by a generalization of this method, one can construct ho- 
mogeneous bodies which are "gravitationally equivalent" , or "gravitationally 
similar" to any given body, or collection of bodies. In Section 4, we build our 
intuition for the resulting shapes of such monopoles by studying, to lowest 
order, some simple examples, in particular the thin cylinder. In Section 5, 
we discuss some of the mathematical properties of the construction, present 
a detailed analysis of the perturbations of a monopole, and outline the pro- 
gram of a rigourous mathematical proof of the existence of such monopoles. 
Physical considerations on the growth of monopoles, including a discussion 

12 As discussed below, our basic method allows generalizations and combinations which 
lead to monopoles of very general shapes, having, for instance, as many holes as desired. 
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of different practical methods for realizing them, are presented in Section 6. 



2 Definitions of monopoles and inexistence of 
"solid" aspherical homogeneous monopoles. 

2.1 Equivalent definitions of monopoles. 

Let p (x) denote a volumic mass density in three-dimensional space. For 
simplicity, we start by assuming that p (x) has a compact support. The 
gravitational potential generated by p (x) reads (with our sign convention, 
and setting Newton's constant to unity) 

U p (x) = j**> 0^ . (1) 

The interaction energy of this mass distribution p (x) in an external gravi- 
tational potential [/ ex t(x) is 

E p = - |d 3 x p(x) £/ cxt (x) . (2) 

Both Up and E p can be skeletonized by introducing the multipole moments 
of the distribution p. They can be described by symmetric trace- free (STF) 
tensors of order I > , 

Qig^M = SIT,;..., /' r/ x p(x) , (3) 

where STF^ i£ denotes a symmetric trace-free projection over the indices 
i\ . . . iif\ and xq the coordinates of some origin O in space. An equivalent 
description uses the spherical harmonics and reads, in the three dimensional 
space: 

Q%o = Jd 3 * p(x) r* YM?) , (4) 

where x l — x l = (r sin 8 cos<p, r sin# siny?, r cos#). 

13 E.g., STFy X' l Xi = X l X? - ±5^ X s X s , 
STF^-fe X l X^X k = X l X^X k - ^ X S X S {S l ^X k + S^ k X l + 5 ki Xi) , etc., in N space 
dimensions. 
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The zeroth order multipole moment (or monopole moment) is simply the 
total mass of the distribution, M = J d 3 x p (x) . We shall always fix the 
origin O at the center of mass of the distribution p, i.e. at the point such 
that the dipole moment Q l p ^ Q = J d 3 x p (x) (x l — x l Q ) vanishes identically. 

A (bounded) gravitational monopole can then be equivalently defined as 
a distribution p satisfying either: 

(i) the external potential generated by p is spherically symmetric: 

UpOx) = | r for x in the exterior 14 of the support of p ; (5) 

|x-x | 

or (ii) the interaction energy of p (x) in any external gravitational potential 
(x) is the same as for a point mass located at the center of mass xq: 

E p = — M f/ cxt (x ) ; (6) 

or (iii) all multipole moments (of order i > 1) vanish identically: 

q = for |>1. (7) 

The equivalence between (i) and (ii) follows from considering the exter- 
nal distribution p ext generating the gravitational potential U cxt . The corre- 
sponding interaction energy is E p [U ext ] = — J d 3 xp(x') p ex t(x) |x — x'| _1 = 
— /<i 3 x U p (x) Pext(x) , to be compared with — M U ex t( x o) = 
— J d 3 x | X ^ Q | Pcxt ( x ) • Demanding the equality (|5[) of these two expressions 
for any p ex t(x) is equivalent to Equation (|5|). 

The equivalence of either (i) or (ii) with (iii) follows from expanding 
|x — x'| _1 in Equation (|I|) in powers of x' — xo , which ultimately yields, for 
x in the exterior of supp p : 

oo (_y I 

u pW = Eir Q ^o d ^-u T——i > (8) 



£=0 



and from expanding C/ ext (x) in Equation (^) in powers of x — x G , which 
yields: 

OO 1 

E p = - E 7T Q%£ d n..M ^t(xo) • (9) 

14 If the support of the bounded distribution p divides space in several regions, we mean 
by "exterior" the outside region, connected to infinity. 
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The properties (i) - (iii) are also equivalent to: 

(iv): the spatial Fourier transform p (k) = / c? 3 x e _ * kx p (x) (choosing 
here O as the origin) may be expanded (around k = 0) as 

p(k) = p(0) + k 2 /(k), (10) 

as we shall discuss later (cf. subsection 4.2). 

Essentially all our results straightforwardly extend to any number of space 
dimensions N > 2 . It suffices to replace the basic Newtonian potential r _1 
by (N - 2)- 1 r"^" 2 ) when N 2 (and by lnl/r for N = 2), and the 
coefficient 4 7r in the Poisson equation 

At/ P (x) = -4 7rp(x) (11) 

by the surface of the unit sphere in R N , ojn = 2 tt n ^ 2 / T(N/2) . To simplify 
the discussion we shall phrase our results by assuming N = 3. 

An homogeneous monopole is defined by a density p uniform within some 
domainF] Q of R 3 . For simplicity, we assume that Q is bounded and admits 
a sufficiently regular boundary S = dQ (assumptions that we shall relax 
later) . 

It is useful to introduce also the concepts of gravitationally equivalent and 
gravitationally similar distributions. Two density distributions p\ (x) and 
p2 (x) will be said "gravitationally equivalent" if they generate the same 
external gravitational potential, i.e. if all their multipole moments coincide: 

QelX 1 = Q£X e for * > • (12) 

(the choice of origin O being arbitrary). The distributions will be said "grav- 
itationally similar" if the external gravity potentials they generate are pro- 
portional. This is equivalent to requiring that their "reduced" multipole 
moments coincide: 

/ M 1 = I M 2 . (13) 

This latter concept is invariant under a (constant) rescaling of the densities 
pi and p2 ■ In this nomenclature, a monopole can be defined as a distribution 
which is gravitationally equivalent, or similar, to a point mass. 

15 By "domain" we mean a connected open subset of M. N . 
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2.2 Inexistence of aspherical solid homogeneous mono- 
poles. 

Let us prove that if Q is a (connected) homogeneous monopole, and if the 
boundary £ = dQ of the domain Q is connected, then Q must be a solid 
sphere. To do so we study the action, on the potential U p (x) generated by 
the monopole Q, of the generators of infinitesimal rotations^ around O, 

I = (x-xo) x . (14) 
ox 

The rotational symmetry of the Laplacian implies that the (Lie) derivative 
of U p (x), say u = L U p , satisfies everywhere, in the sense of distributions, 

A u = - 4vr L p , (15) 

in which the source L p is a "single layer" surface density distributionF]. 

Let us assume that Q defines a homogeneous monopole with a con- 
nected boundary £ = d Q . This implies that £ divides space into two 
regions: the domain Q ("interior of £") and the complement of Q ("ex- 
terior of £"). Each component u % = L l U p of the infinitesimal variation 
of the potential appears itself as a potential, generated by the single layer 
surface density L l p = o % 5 s , distributed on the boundary E . We shall now 
prove that this surface density Uj must in fact vanish identically. Indeed 
the potential u l it generates vanishes everywhere in the exterior of E (since 
Up(x) is rotation-invariant), and is harmonic everywhere inside it. On the 
surface E itself u 1 must be continuous, while its normal derivative might 
a priori have a discontinuity, equal to — 4 n times the surface density a 1 . 

16 If the number of space dimensions is not N = 3, we use 

the subsequent argument remaining the same. 

17 To be explicit, we have 

L p = <7(5 S = - (x - x ) x n p 6-s 

We have used the fact that the gradient of the characteristic function of the domain il is 
— n S £ , n being the outgoing normal to the boundary E = d fl . Ss denotes the standard 
uniform (5-function surface distribution on S (such that J d 3 x /(x) 5% = / s dY. [/(x)]s, 
where dE is the area element on S). 
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The vanishing of u % in the exterior and its continuity across E implies that it 
vanishes on S. It must then vanish everywhere inside E (by the uniqueness 
of the solution of the Dirichlet problem). Since it also vanishes outside, its 
derivative cannot have a discontinuity on the surface £, which implies that 
the density a % distributed on £ has to vanish. The resulting equation, 

Lp = , (16) 

expresses that the monopole Q should be spherically symmetric! 15 ! It must 
be the solid sphere limited by the sphere £ of center O. 

We have not restricted in advance Q to be simply connected. Our result 
therefore excludes, for example, "solid" tori as possible homogeneous mono- 
poles. 

The crucial assumption is that the connected domain Q have a connected 
boundary E (i.e. that we are dealing with a "solid" monopole without any 
internal cavity). By contrast, in the case where the boundary of Q is made of 
two disconnected surfaces, one inside the other (i.e. d Q = S cx tU T, iat , as for a 
"thickened" topological sphere or topological torus, i.e. the volumes between 
two nested surfaces having the topology of S 2 or T 2 ) , or more than twcQ, we 
have no way to conclude that the monopole must be spherically symmetric. 
Indeed we shall see that, in such cases, one can construct infinitely many 
aspherical homogeneous monopoles. 

Our analysis can easily be further extended to aspherical monopoles 
Q = U Qi that would be constituted of several disconnected parts (with 
boundaries £; = dQi). If the exterior boundary E cxt / of one of them does 
not enclose the center O inside it, and is reachable from infinity, it follows 
from Gauss' theorem that the total mass included within it must vanish. As 
a result, the exterior boundary of a gravitational monopole must always be 
connected. For example a collection of solid objects (e.g. solid tori, even 
intertwined) cannot constitute a gravitational monopole. If a monopole is 
constituted of several components all but one should be located within the 
interior of one of them, as in the case of Russian dolls, for example. 

18 Technically, the formula of the previous footnote shows that L p = implies that 
every normal n to £ should pass through O, i.e. that the connected surface E is a 
sphere of center O, and the monopole f2 the corresponding solid sphere. 

19 This is the case of a "Swiss-cheese-like" domain with a disconnected boundary 

s = ar!-E oxt u{E in t I }. 
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3 Constructing aspherical (hollow) mono- 
poles. 

3.1 Thin aspherical monopoles. 

The starting point of our method is a simple consequence of a well-known 
property of electrostatics (see e.g. Ref. fl8fl ). Let Q be a domain with a 



connected boundary £ . Let us consider £ as a grounded conducting surface, 
and introduce a unit (positive) electric charge at some given point O (with 
coordinates x l Q ) within Q . The charge + 1 at O induces on the surrounding 
grounded conducting surface £ a (negative) charge density, say — a, which 
has the well-known property of generating a potential which exactly screens 
the l/|x — xo| potential created by the charge +1, outside £. (This is 
mathematically evident from the uniqueness of the solution of the Dirichlet 
problem A U = , in the complement of Q .) By reversing the sign of 
the surface charge density, we conclude that the positive surface density 
+ a on £ generates, outside this surface, the exactly spherical potential 
+ l/|x-x | . 

In mathematical terms, the potential generated by the charge + 1 in- 
troduced at O is the Dirichlet Green function of the domain Q, G(x, x Q ), 
solution of 

A, G(x jXo ) = - 4tt 5(x-x ) , (17) 

which vanishes when x G £ = d Q . Let us define the following surface 
density on £ 

C7°(x) = -^«9 n G(x,x ), (18) 

where d n denotes the outgoing normal derivative at x 6 £ . The density 
cr° is everywhere positiveF] and integrates to + 1 : 



L 



d£ o-g = - — / rf£ d n G = - — [ rf 3 x A X G = + 1 . (19) 
Air js 4:% Jn 



20 a® is positive even if has a non trivial topology, e.g. that of a multi-handled solid 
torus. Indeed, G(x) being positive (and ~ l/|x — xo| ) near xo , zero on the boundary 
of f2 , and harmonic within Q — {0} cannot (by the maximum-minimum principle) reach 
negative values within Q . Therefore G(x) decreases from positive to zero values when 
crossing E . 
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Applying Green's identity 



/ d 3 x (u A v - v A u) = [ d £ {u d n v - v d n u) (20) 

with it equal to any external potential f/ ext (satisfying A U cxt = within 
HUE), and u(x) = G(x, xo) , yields 

^ext(xo) = / s rfS <7°(x) C/ ext ( X ) . (21) 



Equation (]2lf) expresses precisely the property (y) of a monopole: the 
surface layer with density <t^ interacts with any external potential as if it 
were a unit point mass located at O. (The more "active" definition @ is 
obtained by taking f/ ext (x) = |x — xp| _1 for some external point P.) In 
the mathematical literature, the surface distribution a® c?£ is called the 
"harmonic measure" of £ with respect to the origin O. Equation ([H]) says 
that the value, at a given point O G fi, of any harmonic function can be 
written as a weighted average of the values taken by this harmonic function 
on the boundary of Q . In other words, a® (x), if it is known for all points 
O G Q, gives, via Eq. (pl|), the explicit (unique) solution of the Dirichlet 
problem: determining a harmonic function in Q from its values on £ = dQ. 

In words, the consideration of the surface density a® solves the problem 
of constructing infinitely thin aspherical monopoles. Note that £ does not 
need to have the topology of the sphere. It can have the topology of a torus, 
or of a more complicated surface with many handles. As long as £ is the 
connected boundary of a (sufficiently regular^]) connected open set Q , one 
can construct a thin monopole having the shape of £ . It is surprising, but 
true, that one can choose arbitrarily a complicated multi-handled shape £ , 
and an arbitrary origin O enclosed within £ , and construct, by depositing 
a positive surface layer oc a® on £ , an object generating a potential which 
is exactly proportional to |x — x | _1 outside £ . At this stage, we assume 
that Q is bounded, but we shall later consider the limiting case where Q 
becomes unbounded in a way which ensures that the mass distribution of the 
monopole falls off exponentially fast at infinity. 

21 We refer to the literature on potential theory, e.g. Ref. jl9|], for the characterization of 
the domains admitting a Green function and a harmonic measure. We shall briefly discuss 
later what happens when S is not smooth. 
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3.2 Thick aspherical monopoles. 

To construct "thick" aspherical monopoles, it suffices to build them by grow- 
ing successive layers having a thickness proportional to a® . More precisely, 
let us start with some initial closed surface E which is the (connected) 
boundary of a domain Q and let us choose, once and for all, some origin 
O within Q . It is convenient to use as parameter t measuring the contin- 
uous growing of a thick monopole around X the algebraic volume of the 
constructed monopole, counted positive for monopoles obtained by adding 
layers to the outside of X , and negative when growing E towards the inside. 
At each "time" t the monopole consists of the volume Q t delimited by two 
boundaries: a fixed boundary X and a moving one, say T, t . 

We can (locally) represent the moving boundary S t by equations of the 
form 

x G E t & x = X(t,£) , (22) 

where £ = (£\£ 2 ) denote two (curvilinear) coordinates on the fixed bound- 
ary X . In mathematical terms, X(£, •) is an embedding of X into R 3 . We 
define this embedding (as can always be done) so that points with fixed co- 
ordinates £ propagate orthogonally to the instantaneous X t . The condition 
of adding successive infinitesimal layers of constant density and thickness 
proportional to the harmonic measure a° t reads 

^ X (t, = *° (X(t, £)) n St (X(t, £)) , (23) 

where ri£ t denotes the outward normal to S t . The differential version of 
Equation (|23|), 

rfX = <t£ ( n St dt , (24) 

together with the fact that the surface integral / dY, t cr2 is + 1 (Equa- 
tion (|19D), makes it clear that the infinitesimal volume between T, t and 
is indeed dt , and that each layer d X adds an infinitesimal homoge- 
neous monopole dQ t with fixed origin O. (The constant density p of the 
monopole is set to unity for simplicity.) 

This can also be verified by a simple calculation. Namely, Equations ( p3|) 
or ( pi"D are equivalent to requiring that the t-derivative of the integral of any 
function w(x) over the volume Q t contained between S and S t be equal to 
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the surface integral / d'Et a® tt(X(£, £)). If w(x) is a harmonic function, 
say w(x) = £/ext(x) , the latter surface integral is, according to (0), equal to 
the value of U ext at the origin O. Integrating this result over t we get 

\t\ f/ 0xt (x o ) = / d 3 x f/ cxt (x) , (25) 

which expresses the fact that the volume fl t contained between E and S t 
is a homogeneous monopole, of volume \t\. 

This construction exhausts the possible homogeneous monopoles, only in 
the simple case where: (i) they belong to a continuous^ family of solutions 
with variable mass in which only one of the boundaries of the monopole is 
allowed to move, and (ii) the moving boundary S Q (where a denotes a set 
of parameters) coincides with the fixed one £ for some value of a, for which 
the total mass (i.e. volume 1 1 (a) | ) vanishes. (We assume here, as above, 
that the two boundaries £ and Sck are connected.) 

Indeed, let us consider a (a priori multi-parameter) family of monopoles 
filling the volume Qct contained between the fixed boundary £ and the 
moving one Sq . Let us try to perturb this monopole, at fixed volume t 
(and therefore fixed mass), keeping the boundary £ fixed. The perturba- 
tion of the other boundary Sq is defined by its orthogonal displacement 
5h = (X(a + 5et, £) — X(a, £)) . ns Q . This 5h defines a single layer sur- 
face density distributed on Eq . Since it generates a vanishing potential 
outside this (connected) surface, it must vanish identically. Therefore, the 
only possible variations of S a are those in which the volume (i.e. the mass) 
changes. For these, the uniqueness of the solution of Eq. (pl|), considered 
as an equation for the harmonic measure a, implies that 5h is necessarily 
equal to er^ 8t(at). This shows that, when one boundary S is fixed, the 
other one Sq; must depend only on the single volume parameter t = t(ot) , 
and evolve according to Eq. (]2"5|). Under the further assumption that Eq/ 
coincides with £ when t (ex.) = , we recover the "thickening" process of E 
defined above. 

22 This basic construction method may easily be extended to take into account dis- 
continuous changes, such as, for example, the inclusion of additional cavities during the 
growth of the monopole. There could exist, furthermore, isolated solutions which cannot 
be continuously deformed. 
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This proof of the uniqueness of the contruction of monopoles breaks down 
if we consider continuous families of solutions in which both the outermost 
boundary Eq^ , and one or several inner boundaries E^ 1 are allowed to 
move. We shall see, in subsection 3.3, that the possibility to compensate 
arbitrary variations of Eq;* by appropriate variations of Eq^ allows us to 
construct more general families of monopoles than the basic ones described 
by the thickening process of Eq. (^) - even if this process is generalized 
to allow for discontinuous changes leading, for instance, to the inclusion of 
additional cavities during the growth process. 

Let us mention in passing that one can trivially extend our method to 
constructing many sorts of inhomogeneous monopoles. For instance, if we 
want to prescribe the volumic density p(x) within the monopole it suffices 
to divide the right-hand sides of Equations (^3|) or (p4|) by p (x) when prop- 
agating E t ( \t\ now denoting the total mass of the monopole fit)- Alter- 
natively, one can prescribe the shape of the monopole, and define, from any 
slicing of the monopole by interpolating E t 's a (positive) density p (x) so 
that the orthogonal distance between E t and E t+di is a®Jp. 

Instead of writing the fundamental propagation law for homogeneous 
monopoles in the differential form ([^f), one could as well write it as an 
"eikonal" type equation. Namely, by eliminating the two surface coordinates 
£ 2 ) from the three embedding equations (|2~2"D , we can write an equation 
determining the position of the surface E 4 in the form 

x G E t if (x) = t . (26) 

In this form, the outward normal is n = V</?/|V</?| , and the infinitesimal 
vectors dx connecting E t and Y^t+dt satisfy dt = V<£> • <ix = |Vy?| n • (ix. 
The propagation law dx = andt + tangent vector becomes |Vy>| er = 1 , 

(Vv.(x)) 2 = K x (x)T 2 . (27) 

Contrary to the usual eikonal ( ( Vy?(x) ) 2 = (n(x)) 2 ) or Hamilton- 
Jacobi equations, Equation (|27| ) is not an ordinary local partial differential 
equation for <p(x) because its right-hand side is a non-local functional of the 
function <£>(■) obtained by solving an elliptic problemQ. Correspondingly, 

23 The meaning of crj? x ( x ) appearing on the right-hand side of Equation (|27|) is the 
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the propagation equation ( p3[ ) is local in "time" but non-local in the two 
spatial variables £ 2 . We shall briefly discuss below, in Section 5, in what 
functional spaces we expect the differential equation ( |2"B| ) to admit unique 
solutions (at least for t varying in some interval). 

Finally, let us mention that, when working in any number of space dimen- 
sions^] N >2, the propagation law (^) (or (|27|)) takes the same form. The 
only difference is that S t denotes a moving hypersurf ace (with N—l internal 
coordinates £ ). The harmonic measure a® is still uniquely defined by fl2~T|), 
or explicitly by a® = — ujjy 1 d n G(x, xo) where, as mentioned in Section 
2, ujn denotes the surface of the unit sphere in R N and where the Green 
function has a pole ~ (N — 2) _1 r~( N ~ 2 ^ when N ^ 2 (or In 1/r in N = 2 ) 
when r = |x — xo| — > . 

We shall now extend the basic construction process of Equations (|23|) or 
( P7| ) to more general situations. 

3.3 Gravitationally-similar bodies and more general 
monopoles. 

We have defined, at the end of subsection 2.1, the concepts of gravitation- 
ally equivalent and gravitationally similar bodies, which have identical (or 
proportional) sequences of multipole moments. Let us show how, by an ex- 
tension of our method, we can construct such bodies. A by-product of this 
extension will lead to more general monopoles than the ones obtained earlier 
in subsection 3.2. 

Our electrostatics considerations of subsection 3.1 admit the following 
wide generalization. Let Q denote as above a domain with a connected 

following: Given the point x and the function ip(-), one determines the surface E x passing 
through x as the set of y £ M. N satisfying ip(y) = <^(x) . Then one determines the value, 
at the point x, of the density of the harmonic measure on S x with respect to the fixed 
origin O. 

It is amusing to note that, in dimension N — l, any bounded object, say any (finite) 
reunion of disjoint finite intervals, constitutes a homogeneous monopole. Indeed, a general 
harmonic function is simply u(x) = a + bx so that the condition (Q) holds simply by 
linearity. Still, if we consider the analog of our general construction, i.e. the growing of 
two intervals located around their fixed center of gravity, one finds that the dynamics of 
their growth keeps (in a simplified form) some of the features of (p3|). 
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boundary E . Let us consider E as a grounded conducting surface, and 
introduce within it an arbitrary (fixed) volumic charge distribution — p\ (Q) , 
where Q denotes a generic point of Q . This charge distribution induces on 
the grounded boundary E a surface charge density obtained by superposing 
all the elementary charge densities induced by point charges at the running 
point Q. If cr^ (x) denotes the surface density ([TJ) taken for a pole O located 
at Q (see also Eq. (|77f) below), the resulting surface density induced on E 
by — Pi (Q) reads: 

<t#(x) = JdV Q Pl (Q) <rg(x) . (28) 

By construction, this surface density generates, in the exterior of E , the same 
potential as the arbitrary volumic distribution p\ contained within E . In 
other words, the surface distribution cr^ 1 (x) is gravitationally equivalent to 
the volumic distribution + pi (Q) . 

To construct an homogeneous 25 volumic distribution p 2 equivalent to pi , 
we can now turn on continuously the distribution p x (by considering t p x , 
with t increasing from to 1), while growing, at the same time, successive 
surfaces E t from the initial surface E , in such a way that the elementary 
volume d Q t of density p 2 enclosed between E t and E t+ ^ is gravitationally 
equivalent to the infinitesimal mass distribution pi dt . This new growing 
process is defined by the following differential equation (which generalizes 
Eqs. ©or©): ' 

dX = — a£ n St dt , (29) 

P2 

where n St denotes the outward^ normal to E t , and the surface distri- 
bution ( p8|) taken for the moving surface E t . The total mass of the corre- 
sponding layer dVt t - which by construction generates at the exterior of E t 
the same potential as the infinitesimal distribution p\ dt - is easily verified 
to be (as it should) 

dM 2 = p 2 I d E t n St . d X 

25 The construction given below trivially extends to inhomogeneous distributions. 

26 To fix ideas we consider here an outward growth process, which has the advantage 
of being stable (as discussed later); we can define as well an inward growth process, by 
changing the sign of the right-hand side of Eq. (^9|) . 
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= f dZ t erg] dt = I dV Q pi(Q) dt = M l dt . (30) 

By ending the growth process at tf — 1 we obtain a thickened version of 
£ which, when filled with density p 2 > is gravitationally equivalent to the 
volumic distribution pi . Distributions gravitationally similar to p\ can be 
obtained by a rescaling of p 2 , or by ending the growth process at tf ^ 1 
(thereby generating a distribution gravitationally equivalent to tjpi ). 

This construction is interesting in several respects. From a theoretical 
point of view, it solves a generalization of the monopole construction prob- 
lem by providing, for any (connected) surface £ and any mass distribution 
Pi located within £ , an homogeneous thickened version of £ which is grav- 
itationally equivalent (or similar) to p\ . (Note that pi, being totally arbi- 
trary, can correspond to a collection of disconnected bodies with arbitrary 
densities.) 

It also leads to the construction of more general monopoles than the ones 



obtained by the basic growth process ( 23j , ^), even when it is generalized by 



introducing discontinuities, with the inclusion of additional cavities during 
the growth. (This generalization can be viewed as the superposition of several 
"basic" monopoles, all of them having the same center of mass O.) 

To show this we start from an homogeneous monopole M. of center O 
(e.g. a solid sphere, or one of the monopoles constructed above). Let us 
dig arbitrarily various holes within this monopole. The mass distribution 
Pi (Q) that we take off corresponds to a subset 7i of M. . This density 
Pi (Q) taken away from the initial monopole is gravitationally equivalent, as 
we have seen, to a new homogeneous body AA4 (filled with the same density 
as .M), obtained by thickening, towards the outside, the outer boundary E of 
our initial monopole M.. Altogether the resulting object {{M. — 7~L) U AM.) 
is a new homogeneous monopole, of the same mass M and center of mass O 
as the initial one. But it now presents a new set of holes TC , corresponding 
to one or several disconnected volumes of arbitrary shapes. 



The construction ( p9|) is also quite interesting from a practical point of 
view. One motivation of this work is to define differential accelerometers 
which are - ideally - totally insensitive to arbitrary external gravity gradients. 
This may be realized by using two gravitationally similar bodies centered at 
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the same point, the accelerations induced by the external gravity gradients, 
proportional to QV^^* / M Pi , being identical. 

Given some inner body with mass distribution pi (Q) (e.g. a solid straight 
cylinder) and an arbitrary surface E enclosing it, we can now define, from 
the growth process (^9]) considered with any value of the density pi and any 
final "time" tf, a one-parameter family of homogeneous outer bodies gravita- 
tionally similar to the inner one. These bodies have the form of homogeneous 
thick shells built on the boundary £ . We shall see later how this construc- 
tion can be extended to the case where E has, for example, a cylindric or 
cylindriclike shape. 

Before discussing in Section 5 some of the mathematical properties of 
the evolution equations ( ]2"5| ) or fl2T3|), we shall build up our intuition on the 
growing of monopoles by solving explicitly, in Section 4, the first step of the 
construction in some simple geometries. 



4 Simple examples of thin monopoles. 

4.1 The plane and the sphere. 

The most trivial example would be the limiting case where the domain Q 
is a half of M 3 , £ = dfl being an infinite plane. The Green function with 
pole at O G Q is simply G(x, x G ) = |x — x | _1 — |x — Xq| _1 , where O is 
the image of O through the "mirror" £ . Equation ( |T5D yields the surface 
distribution on £ , 

°s — 3 ' (31) 

2 71 T 

where r = |x — xq| and a is the orthogonal distance between O and the 
plane E . As seen from the other half-space that is the complement of Q , 
the surface distribution (|3lD generates exactly the potential |x — xo| _1 . 
However, as the asymptotic fall off of a is rather slow it does not make sense 
to compute the multipole moments of the corresponding surface distribution, 
nor to say that we have generated an interesting thin monopole. 

The second most trivial example consists of the case where Q is a ball 
of radius R in K 3 , E = dQ being a sphere. Let C be the center of the ball 
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and O be any point within £7 . The Green function is again given by the 
method of images: 



G(x,x ) = Ix-xoP ~ AIx-XqP where A = R/CO , (32) 

and where O is the inverse of O through the sphere E (i.e. CO = R 2 CO/ 
|CO| 2 ). The corresponding surface distribution is given by the well-known 



Poisson formula (see e.g. Ref. |T9|j ) 



o R 2 - CO 2 1 

= ArrR ^ ' (33) 

which exhibits the same r~ 3 dependence (where r = |x — xo| with x 6 E ) 
as the planar distribution (^T|). (In iV space dimensions the Poisson formula 
is simply obtained from ( |3"3"D by replacing Att — > ujn and r 3 —> r N .) 

We can now consider the surface distribution (|33"D on this sphere (of ge- 
ometrical center C) as defining a thin monopole (of center of mass O). By 
construction, it generates a Newtonian potential which is exactly |x — xq| _1 
outside E , its interaction energy with an external potential is exactly that 
of a unit mass point located at O, and all its multipole moments (of order 
£ > 1) around O vanish identically. 

The lesson we learn from ([y]) and (|33D (at least when one starts with 
a surface E with a simple shape varying only on length scales larger than 
the distance between the pole O and E ) is that the growing of a monopole 
will consist of successive layers on E forming bumps roughly concentrated 
around the points of E which minimize the distance |x — xo| • When grow- 
ing the monopole towards the outside, one expects the bumps to become 
smeared (because their centers move away from O), and the surface E t to 
asymptotically approach a very large sphere of center O. On the contrary, 
when growing E t towards the inside, one expects the bumps to start de- 
veloping into fingers rapidly advancing towards O. We shall further discuss 
later the instabilities associated with the inward evolution and the smoothing 
character of the outward propagation. 

4.2 The thin cylinder. 

A less trivial - and much more interesting - example of a thin monopole is 
obtained by choosing for surface E = d Q a straight cylinder (of radius R) in 
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M 3 . We choose, for simplicity, the pole on its symmetry axis. In cylindrical 
coordinates ( p, ip, z) the pole is at p = z = and the equation of £ is 
p = R. There are (at least) two equivalent ways of finding the corresponding 
Green function. Let us start with one before explaining the other. 

Writing the Green function G?(x) = - + U (p, <p, z) , where r = y/ p 2 + z' 2 , 
the function U (manifestly ^-independent) must be harmonic within O , 

d 2 U IdU d 2 U 
op 2 pop a z 2 



and equal, when p = R, to —1/r = — l/\/R 2 + z 2 . Taking a Fourier 
transform with respect to z, U (p, z) = dk (p) e tkz , we find from 
fl3"4|) that Uk{p) is proportional to a modified Bessel function of order zero: 
Uk(p) = Cfc Io(kp) ■ The coefficient is obtained by writing that Uk(R) 
is the inverse Fourier transform of U (R, z) = —1/ \J R? + z 2 . The lat- 
ter transform is given by a modified Bessel function of the second kind: 
/ °° dz coskz/Vz 2 + R 2 = K (\k\R). Finally, we get 

1 , rrf. 1 1 f + °° ,u M\k\R) TtU^JO* 



z) = 7WTT 2 + f/(p ' z) = 7WTT 2 ~ « L dk -htm /o(M 6 

(35) 

The corresponding surface density (|18D can be written (when writing also 
the (p 2 + z 2 )~ l l 2 term as a Fourier integral with respect to z, and using the 
Wronskian k [Iq(k) K (k) — Io(n) K' q (k)] = 1 ) as 

where (using the dimensionless variables ( = z/R and k = fci?) 

/+oo p* K C 
S 77^ • (37) 

A different method of getting G and a is to take advantage of the 
vanishing of G ( p, z) at p = R (for any z) to expand the Green function 
G ( p, z) as a Fourier-Bessel series 

2 00 

G = ~5 H a e (z) J (k e p) , (38) 

^ 1=1 
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where kgR = ~ 2.405, 5.520, 8.654, . . . are the successive zeroes of the 
Bessel function Jo . Inserting this expansion into Poisson's equation AG = 
— 4 ii 5 (i.e. Eq. ( |3"4"D completed with the source term — 4 ix S(z) S±(p) ), one 
gets: (i) the information that ae(z) = Aie~ k ^ , and (ii), from the strength 
of the delta function, the values of the coefficients At . (More details on this, 
and on the properties of the resulting distribution function F , will be given 
elsewhere.) Finally, we get 

2 00 1 

* l| ^^ , "" l,lM ' (39) 

in which enters the Bessel function of order one, J\ (x) = — J' Q (x) , evaluated 
at the zeroes of Jq{x) . From ( p9|) one deduces that the surface density (|I~8|) 
can be written as in Equation (|3~6"D, with 

OO 1 

HO = E T7-T ^ ICI • (40) 

The identity of the different-looking results ([J7|) and (^) is easily verified 
by folding (when, say, ( > ) the k contour of integration in (|37| ) (after 
closing it by an infinite half-circle in the upper complex k plane) around 
the upper imaginary axis k — % x , thereby picking up an infinite series of 
contributions due to residues at the simple poles of 1/I (ix) = 1/ J (x) . The 
Fourier-integral representation (|37D is valid for all values of ( and makes 
it easy to compute numerically the shape of the function F((). The series 
( [40"D is absolutely convergent only when C ^ [the slow decrease of Ji(je) ~ 
(— y +1 i^ 1 / 2 as £ — > oo causes the series (|40"1) to diverge when ( = 0], but 
gives a useful representation of F(() when \(\ ^ 1 because it captures well 
its asymptotic fall off. In particular, the leading behaviour for \(\ > 1 is given 
by the first term in PO]), coming from the first zero of Jo(x) : j\ ~ 2.40483 
(and Ji(ji) — 0.51915). Therefore, for \z\ ^ R the surface density (|36|) is 
approximately given by 

~ e - 2.405 M ( } 

The function F(Q (which integrates to unity) is represented graphically in 
Figure 1. 
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By construction the surface density (|36| ) laid on an infinite straight cylin- 
der of radius R generates a l/|x — xo| potential everywhere outside the 
cylinder. This monopole is of infinite extent, but, thanks to the exponential 
fall off ( flip of the density it makes sense to consider its multipole moments 
of arbitrarily high order, all the integrals / d'E o STF(x n . . . x H ) being con- 
vergent. It is clear (by considering, for instance, the infinite cylinder as the 
limit of a finite cylinder closed by two caps) that our construction guarantees 
that all the multipole moments of the distribution ([36]) vanish identically (for 
£>1). 

This can be verified by explicit computations. A very simple way to 
proceed (which extends to all distributions which are bounded, or have an 
exponential fall off) is to consider the spatial Fourier transform of the density 
distribution: p (k) = / d 3 x e _tk ' x p (x) . This is an analytic function of k 
which can be analytically continued to (sufficiently small) complex values of 
k. We define the "inertia moment" of order i of p(x), say J* 1 as the 
integral on the right-hand side of Eq. (E3) without the STF projection. These 
inertia moments are the Taylor coefficients of the (convergent) expansion of 
p (k) in powers of k. Each inertia moment P 1 - 1 * is equal to its trace-free 
projection Q n — ie plus some "trace terms" containing at least one Kronecker 
fii a % _ jt j s ea sily seen that a monopole (Q n ~ H = 0) can be characterized 
by requiring that all terms of the Taylor expansion of p (k) contain a factor 
k 2 = 5ij k l W . In other words, a monopole is characterized by requiring that 
the Taylor expansion of p (k) reduce to p (0) when the complexified k is 
only restricted to satisfy k 2 = . This general characteristic property of a 
monopole can also be proven by taking the Fourier transform of the Poisson 
equation satisfied by the difference potential V(x) = U(x) — — , where U(x) 
denotes the potential generated by the monopole (centered at the origin) and 
M its mass. Indeed, A X (U-^-) = -An (p (x) - M5(x)) , which implies 
p (k) = p (0) + 7— k 2 V"(k) , where V(k) denotes the Fourier transform of the 
compact- support function V(x). 

For the case at hand, the explicit calculation of p(k) = 
J dH er(x) e _lk ' x , with k = (/icos0, /isin0, k) , is very easy using (|3~7l) 
and yields p (k) = Jo (h)/Io (k) . The fact that it reduces to p (0) = 1 when 
h = ik (i.e. k 2 = h 2 + k 2 = ) establishes directly the vanishing of all the 
multipole moments. The same verification can be done from inserting (|40|) 
in (^) at the price of more Bessology. 
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The important lesson we learn from this last example, by truncating a at 
a finite \z\ > R , is the possibility to construct exponentially- accurate approx- 
imations to exact monopoles having the shape of finite cylinders, opened at 
both ends, and thickened in the appropriate way in the middle. From (ffl]) we 
see that the accuracy with which the multipole moments of such a truncated 
thin cylinder vanish is of order e~^ z " R . More generally, one can get such 
exponential decreases of a at infinity, whenever, when some longitudinal di- 
mension (playing the role of z) tends to infinity, the transverse dimensions 
remain boundedF]. 

Such bodies could be very good candidates for being the "external cylin- 
der" of a differential accelerometer of the STEP type. They are exponentially 
insensitive^ to external gravity gradients, and allow easily for an "internal 
cylinder" (together with its suspension and sensing system) to be placed 
inside them. 

Alternatively we can use the generalized construction of subsection 3.3 
to define, for a given internal straight cylinder of finite length - with non- 
vanishing multipole moments - a corresponding "external cylinder" having 
exactly the same reduced multipole moments Q nl2 -^ j M . In the thin cylin- 
der approximation, the appropriate surface distribution a is obtained by 
generalizing the solution (|35| - ^0D to the case where Eq. ( pi[) is replaced by 



A G p (x) = — 4 7r p (x) , where p (x) denotes the mass distribution of the 
internal cylinder. The exponential fall-off of o will be given by the same 
exponential factor as in Eq. ( f4~0l) (with, in general, a different prefactor). 
We shall come back, at the end, to the possible use of this construction for 
Equivalence Principle experiments. 



27 The asymptotic behavior of a is then primarily determined by studying the eigenval- 
ues of the Schrodingerlike equation — At ra nsvcrsc ip — E ip within the interior of £ , with 
Dirichlct boundary conditions. The crucial criterion for an exponential decrease of a is 
the existence of a positive lower bound on E (typically proportional to the inverse of the 
square of the transverse dimensions). 

28 The suppression factor is roughly proportional to exp(— ji L/2R), where L is the 
total length of the cylinder. 
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5 Mathematical aspects of the construction 
of aspherical monopoles. 

5.1 The mathematical setting. 

We have shown in Section 3 that the problem of constructing aspherical 
monopoles with a given center of mass0 amounts to solving a differential 
equation in the functional space describing closed^] surfaces £ . It seems 
intuitively obvious that, if one starts from a sufficiently smooth surface So 



at "time" t — , the differential propagation law (|24|) , <iX = a^ t ri£ t dt , will 
generate, at least for a finite interval of values of the parameter t around zero, 
a continuously evolving family of neighbouring surfaces S t . The problem 
of characterizing the class of good initial data S which admit such an 
evolution is, actually, highly non trivial. From the considerations we shall 
now present, it is clear that (if we demand evolution for both signs of t 
around zero) we need to work with analytic (C w ) surfaces embedded in K . 
We have, however, not precisely determined what class of norms one must 
use to guarantee that S is "sufficiently analytic". We shall see also that 
the problem of outward propagation (t > 0) is very different from the inward 
propagation one, and allows one to start with much less regular (say, only 
continuous) initial dataQ. 

A first clear hint that one needs to work with analytic data is the fact 
that, in two dimensions, the propagation law (|23"D is reminiscent of the 
Cauchy-Riemann equations for the mapping (t, £) — > X(t, £) = (X, Y) , 
dX/dt = <9F/<9£, dY/dt = —dX/d^, which can be written as 
dX(t,£)/dt = |<9X/d£| n(X). Clearly the local functional |0X/d£| is 
a much simplified version of the non-local functional o"s(X) appearing in 
( |23|) which has, as we shall see, a worse sensitivity on the regularity of the 
function X(f) than \dX/d£\. 

29 In this section, the origin O used in the construction of Section 3 is fixed from the 
beginning. We henceforth drop the superscript O on a§ . 

30 By "closed" , we mean more precisely that £ is the connected boundary of a connected, 
open subset of M. N . 

31 The fact that the evolution is stable outwards and unstable inwards has been known 
for some time in a very similar two-dimensional problem (injection of fluid in a Hele-Shaw 
cell; see e.g. Ref.[p0|). 
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Having recognized the need to work within a (real) analytic setting, we 
can abstractly formulate our problem as a functional equation (for u belong- 
ing to some functional space) of the form 



du(t) 
dt 



F[u{t) 



(42) 



which is the ultimate generalization of the Cauchy-Kovaleskaya problem. 
Nirenberg [21] and Nishida P2 have proven a general theorem about the ex- 
istence and uniqueness of solutions to the general abstract non-linear Cauchy- 
Kovalevskaya problem (0). The crucial point of their work was to formu- 
late the abstract analog of the condition (in the usual Cauchy-Kovalevskaya 
problem) that the functional T [ u ] on the right-hand side of Equation (fi2]) 
involve at most first-order (spatial) partial derivatives of u and be linear in 
them. Their formulation assumes that u varies in a "scale of Banach spaces" 
S = { B s ; < s < 1 } such that, B s C B s r, when s' < s , with the norms 
satisfying || • || s / < || • || s . Then the crucial "quasi-linearity" condition (ex- 
pressing the analog of the linearity in first-order derivatives) is that T [ u ] 
must satisfy a Lipschitz-type inequality: 



F[u] -F[v\ \\ s . < 



C 



u — v 



for s' < s 



(43) 



Under ([|3|), and some other mild assumptions, the solution of ( f42"D exists 
and is unique (for sufficiently small t). The norms || || s appropriate to our 
problem are the ones standard when dealing with analytic functions (e.g. 

oo oo 
n=0 s n=0 

embedding functions £ G £ — > X G M. N 



or slight adaptations of these to deal with the 



Actually, the functional of X(£) appearing on the right-hand side of 



the propagation equation (g3|) cannot satisfy the condition ( [430 because it 
contains the product of two "dangerous" functionals: the normal (which 
depends on the first derivatives <9X/<9£) and the harmonic measure as . In 
order to apply the abstract theorem of Refs. p2| , one should transform 
the original evolution equation (^) into an equivalent "quasilinear" system 
that we shall introduce later, in subsection 5.6. But, first, to understand the 
subtleties of Equation (|23|) let us discuss the linearization of our evolution 
problem, by studying the small perturbations of a monopole. 
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5.2 A linear propagation equation for the perturba- 
tions of a monopole. 

In this subsection we discuss how a monopole, built from an "initial" surface 
S , reacts to a small deformation of this surface S . The deformation 
of S is parametrized, at first order, by a "height function" h(t — 0, £) 
defined on S . This deformation, which extends to the whole monopole, 
i.e. to the whole family of surfaces £j, is described - still at first order - 
by a function h ( t, £) which gives the local orthogonal displacement of the 
intermediate surface S t . The quantity h(t,£) obeys a linear propagation 
equation, in terms of the "time" evolution variable t (i.e. the algebraic 
volume of the monopole Q t included between the two surfaces S and T, t ). 
We shall establish below, by different methods, this evolution equation of 
h, and discuss, in subsection 5.3, some general features of its solutions. In 
subsections 5.4 and 5.5 we shall explicitly solve it in the case of the outward 
propagation, and discuss, as an illustrative example, the simple case for which 
the initial (unperturbed) monopole has a spherical symmetry of center O. 

The function h ( t, £) will be shown to satisfy a linear (real "Schrodinger- 
type" or "heat-type" ) propagation equation, 

^ h(t) = - H(t) h(t) , (44) 

in which the (self-adjoint) operator Tt(t) is the one defining the variation 
of the harmonic measure <j£ t under the perturbation h(t) , according to the 
equation 

Sax t = - H(t) h(t) . (45) 
This "Hamiltonian" H is explicitly given by 

H : h -> Hh = A/" St (<7 Et h) + ka^h , (46) 

in which k denotes the mean curvature of the surface S t (i.e. + , in 
three dimensions), and As t is the "Neumann operator" which associates, 
to every continuous function (u) on T, t , the value (d n u) of the normal 
derivative of its harmonic extension (within the inside of £ t ). (This propa- 
gation equation of the deformation h , incidentally, must admit as a solution 
<7£ t itself, since h(t) = a^ t 5t corresponds to an obvious and rather trivial 
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deformation of the monopole associated with an infinitesimal shift 5t in the 
variable t .) The essential positivity of this Hamiltonian (meaning that it can 
be made positive by adding a constant) will imply that the outward (t > 0) 
propagation of the perturbation h will be well-defined and smoothing for all 
t , as we shall discuss in subsections 5.3 and 5.4. We shall also solve explicitly 
this propagation equation ([0]), in the case of the outward propagation. 

Let us now assume, after this general presentation, that we know some 
exact solution of the evolution problem ([23|), defining a monopole. (At this 
stage, the only exact solution we can write down easily corresponds to the 
trivial case where So is a sphere and the origin O is chosen at its center, the 
T t 's constituting concentric spheres parametrized by the algebraic volume 
between So and T t .) We consider a small perturbation of this solution. If 
X(£, £) denotes the embedding associated with the unperturbed solution St , 
let X'(i, £) = X(£, £) + eXi (£,£), with e some infinitesimal parameter^, 
denote the embedding associated with some neighbouring solution T[ . Note 
that we consider here variations of end surfaces, (S , T t ) — > (S' , T' t ), keeping 
fixed the volume t spanned between the surfaces. 

Modulo a (t-dependent) reparametrization of T' t , the perturbation of T t 
into T' t can be described by a simple scalar quantity: the orthogonal distance 
between T t and T' t (i.e. the local "height" of T' t with respect to T t ), say 
(after division by e) 

h(t,£) = n(t,£) ■ Xi(U) . (47) 

This scalar quantity also appears naturally if we perturb the representa- 
tion (PB|) instead of the embedding representation fl2"2"|). More precisely, if 
<£>(x) = t is the equation of the unperturbed family of surfaces T t , and 
<£>'(x') = <£>(x') + e v?i(x') = t that of the perturbed family T' t , one finds 
h = - (fx/] V(p\ = - cr(x) ^i(x) . 

An evolution equation for h can be obtained by perturbing (and project- 
ing) the vectorial equation (^3|), or, equivalently, by perturbing the (scalar) 
eikonal equation Q2"Tp. To describe the result we must first discuss some math- 
ematical objects entering the change of the harmonic measure <te under an 
infinitesimal variation of the surface S . 

32 Henceforth, we work to first order in e, disregarding all terms of order e 2 , ... . 



29 



Let us first introduce the "Neumann operator" As associated to S = 
dfl, where Q is some domain in R N . The functions which are harmonic 
within Q (and continuous on its closure Cl = Q U E ) are uniquely determined 
by their values on the boundary E = <9fi ( "Dirichlet problem" ) . If £ denote 
some coordinates on E , any sufficiently regular (say continuous) function 
u (£) on E defines a unique harmonic function u (x) within (Aw = 
0), called the "harmonic extension" of u(£) . Under some mild regularity 
conditions for u (£) and E , the normal derivative of w (x) on E exists. 
Evidently, [<9 n w] (£) is a linear functional of u (£) . This linear functional is 
called the Neumann operator: 

[d n u}v = As « • (48) 

This operator As acts within the space of (sufficiently regular) functions 
on E . It is a positive, elliptic pseudodifferential operator of order 1, which 
means essentially that its kernel A"(£, £') is singular as r] = £ — — > in a 
manner such that its Fourier transform with respect to rj is of order |k| as 
|k| — > oo (k denoting the momentum associated with 77). It is also known 
that the leading short-distance singularity of Af is given by y '— A^ where 
As is the intrinsic Laplacian on E , as defined by the inner metric induced 
from the ambient euclidean metric on R N (see e.g. Ref. [p3fl ). 

By using the second Green's identity 

/ d E («a„ti-w4M) = / d N x (fi A v - v A fi) = , (49) 
js 

one sees that the Neumann operator is self- adjoint with respect to the scalar 
product (u, v) = J-^dT, uv between functions on E: (u, Afv) = (A/"w, v) . 
Moreover, J\f is a positive operator, with respect to this scalar product, 
because of the first Green's identity 

/ dE u d n u = f d N * V{uVu) = f/x [(V«) 2 + «A«], (50) 

applied to harmonic functions: 

(u, jVu) = [ d N x (Vw) 2 > 0. (51) 
Jn 

Note that the positivity of J\f has resulted from using the outward normal 



in its definition (4"8"D 
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The operator A/s determines the modification of the harmonic measure 
erf? under an infinitesimal variation of the surface E . Let us describe, as 
in Equation (fT7|), an infinitesimal variation of E by the orthogonal distance 
between E and E' , i.e. as the map 

(X(£) G E) - (X'(0 = X(0 + e M^n G E') . (52) 

We can compute the change 8 a® (£) = cr°,(^) — er°(£) (evaluated, following 
( |52| ) , at two orthogonally corresponding points on E and E' ) by varying the 
defining relation 

u(O) = JjY> «(X(£)) , (53) 

in which u(pc) is an arbitrary harmonic function. The variation of the left- 
hand side of fl53D vanishes, since one keeps both the origin O and the function 
u (x) fixed. The variation of the right-hand side of ([53]) yields three contri- 
butions. 

The first one is due to the change of the area element d E (corresponding 
to a given coordinate span d N_1 $, ) under the orthogonal deformation ([52]). 
This is described by the trace of the second fundamental form of E . We 
recall that, using some Gaussian coordinates, the metric of R N around E 
can be written as (a, (3, 7 = 1, . . . , N — 1) 

ds 2 = dn 2 + g af3 (n,C) dCdf 3 , (54) 

where n = is the equation of E (so that g a p is the first fundamental form 
of E) and where the curves £ Q = const, are ambient geodesies (i.e. straight 
lines in R N ) orthogonal to E . The second fundamental form (or extrinsic 
curvature tensor) of E is defined as one-half the normal (Lie) derivative of 

9a/3 '■ 

Mn = \ 

For a surface in R 3 , the diagonalization of the mixed tensorf^] k a ^ = g ai fc 7j g 
defines the two principal curvatures k\ = 1/R± and k% = I/R2 ■ From 

33 The mixed form of k can also be defined as the "Weingarten" map, Vv n = fc(V) , 
where V is a vector tangent to E , and [fc(V)] a = k a (i V 13 . 



dg a p{n,C' 
dn 



(55) 



. n=n 
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d E = (det g a p) i / 2 d N - x i it is easy to see that the variation of d E under the 
orthogonal deformation ( j52l) depends only, at first order in e, on the local 
value of h(£) (but not on its derivatives). It is given by 

dZ' = rfE [l + e h(£) k(£)] , (56) 

where 

k(£) = k a(3 ee k a a (57) 

is the trace of the second fundamental form, also called the "mean curvature" . 
In the most relevant case of a surface in E 3 , k — 4- + 4- . 

Coming back to the variation of the right-hand side of Equation (|5*B]), we 
find the following three contributions 

= J^dE [(ehk)u + 5a% u + {eh d n u)\ . (58) 

The third contribution (which comes from the variation of u (X(£) ) caused 



by the change (|52"D of its argument) can be rewritten in terms of A/s u . 
Finally, using the self-adjointness of the Neumann operator, we can rewrite 
( |58|) as a linear form in u (£) . The arbitrariness of u (£) then yields the 



looked for formula for the "normal" variation of a® 

8 = = -e[M s + k}(a<ih) . (59) 

Here, Af(ah) = J <2E(£') A/(£, £') a(g) h(£') is a non-local functional of h , 
while ka h = k($,) o~(£) h(£) depends only on the local value of h at £ . (The 
operator As + k is independent from the choice of the pole O; therefore, by 
linear superposition, the equation ([59]) gives also the infinitesimal variation 
of the weighted harmonic measure a% defined by Eq. (|28|) and entering our 
generalized construction (|29]).) 



We are now (nearly) in measure of deriving an evolution equation for the 
perturbed "height" h(t,£), defined by Equation fl4T|). A direct perturbation 
of the propagation law (|^) yields 

e dX f t '® = o§{X') n'(X') - og(X) n(X) . (60) 
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Projecting this vectorial equation along the normal n yields (to first order 
in e) the following expression for the evolution of the height h — n • Xi 

^ = ^( n " X i) = <ri,(X + eX 1 )-ag(X) + ^.X 1 . (61) 

Because the displacement Xi has both a normal and a tangential compo- 
nent, Xi = nh + Xt , the variation of cr^ appearing on the right-hand side 
is the sum of two contributions: (1) the term fl5"9|) due to the orthogonal dis- 
placement h = n • Xi , and (2) the supplementary contribution + Xy • V a° 
due to the tangential displacement. However, the latter is precisely can- 
celled by the last term on the right-hand side of Equation fl6~ip . Indeed, the 
t-derivative of the normal field n(£, £) is given by 

|^n(t,0 = -V T a° , (62) 
where denotes the tangential gradient (along H t ). 



One way to prove (p2|) is to notice that, by definition of the t and £ 
parametrizations, the metric element of E , ds 2 = <iX 2 , has the form, in 
(t, £) coordinates, 

ds 2 = a 2 {t,C)dt 2 + g a p{t,C) dC df . (63) 

Writing the curvilinear components of the absolute acceleration of the lines 
£ 7 = const, is very easy in the coordinate system ( |63| ) (only one Christoffel 
symbol survives) and yields (d 2 ~K/ds 2 ) a = (dn/ds) a = — d a ojo , which is 
equivalent to 



Altogether, using ( |59|) we get from flSl] ) the simple evolution law 

§- t Ht,£) = -\N+k]{oh). (64) 



Apart from h , all quantities in (|64|) refer to the unperturbed embedding 
x = X(£, £) and its associated family of surfaces E$. The explicit expression 
of the right-hand side of (Bl) may be written as follows: 



" d N ~ l i' Jgtf) Af&g) a(0 h(t,£) - k(i) a(£) h{t,Q , (65) 
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where g = det g a p . The evolution law ( pip applies also to the perturba- 
tion of the generalized growth process ( p9|) leading to the construction of 
gravitationally-equivalent bodies, if we interpret a as the weighted harmonic 
measure cr^ 1 of Eq. (|28"D . 



An alternative way of deriving (|64"D consists of perturbing the (scalar) 
eikonal equation (|27|). Inserting 

<^'(x) = (p(x) + £<^l( x ) an d °E' ( x ) = a ( x ) + £°"l( x ) (66) 

into (^) yields V tp • V = — cr~ 3 o - ! . The perturbed harmonic measure 
is obtained from (^) with the subtlety that the height h is not the expected 
— (fx/l'Vcpl = — a i , but instead h(t,£') = — <x(£') [ (pi(£') — v^i (€) ] 5 the 
last contribution coming from the fact that, in this formulation, the surface 
S x is constrained to pass through the point x when <p — > ip 1 (so that the 
volume parameter £' = <//(x) 7^ t = ip(x.) ). The resulting equation yields back 
( |64|) upon using the following expression for the t-derivative of the harmonic 
measure: 

§- t o{t,t) = -[N + k\{a*) . (67) 

The latter result follows directly by applying the general result ( |59|) to the 
case where S = S 4 and £' = are two neighbouring surfaces within 

the unperturbed family of surfaces satisfying (^J). Note that Equation (|5"7j) 
says that h = a is a particular solution of the general perturbation equation 
(|6"4"1). [This is indeed, as we have already mentioned, a "trivial" solution 
corresponding to considering the t-shifted embedding x = X(t + e, £) as a 
"perturbation" of x = X(t, £) •] 

A third, instructive, way of deriving the perturbed evolution equation (|64]) 
consists of directly varying the fundamental property defining a monopole, 
i.e. 

1 1 1 u (O) = / e^x w(x) , (68) 

where 'S(x) is an arbitrary harmonic function, and fl t the domain (of volume 
1 1 1 ) extending between the initial surface S and E 4 . Keeping fixed the 
origin O, the volume t and the harmonic function w(x) , the change of (|68|) 
under variations of the two end surfaces Sq and S t (both described by a 
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normal displacement as in fl5"2|)) is easily seen to lead to a "conservation" P] 
equation = e (It — Io) where 

It = [ dE t h(t,x) u(x) (69) 

is a surface contribution associated with the variation of the "upper cap". 
Passing as above to the (unperturbed) variables (t, £), one finds easily (using 
the tools given above) that the "time" derivative of I t can be rewritten as 

% = L dE **{Tn + w +k K"o} ( = ")■ < 70 > 

thereby recovering (|64|). 



5.3 Behavior of the solutions of the propagation equa- 
tion for the perturbation of a monopole. 

One can consider Equation (|6~4l) as a Schrodinger equation for the "wave 
function" h(t, £) with respect to the imaginary time t (or, directly, as a 
type of "heat equation"). To do this, it is convenient to endow the space 
of functions on a surface £ with a new scalar product differing from the 
canonical "metric" one [ (u, v) = J^dT, u v] used above. Let us define the 
"harmonic" scalar product (for some given origin O) 

(u,v) = [ d£ a% uv . (71) 

We also define the following operators acting on the space of functions on 
some surface E, 

Ho : h -> H h = AT (a h) , 

(72) 

H k : h -> H k h = (H + ka)h = [Af + k] (ah) . 

The important point is that both Ho and Hk are self-adjoint with respect 
to the new scalar product (|7T1) . Indeed, (u,Hov) = (a u, Af(cr v)) = 

34 By "conservation" we mean t-independence, as in classical mechanics when t denotes 
the time. 
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(J\f(au), av) = (TCqu, v) . Another important point is that the positiv- 
ity of Af with respect to the "metric" scalar product, implies the posi- 
tivity of Ho with respect to the "harmonic" scalar product (]TT|). Indeed, 
(u, H u) = (cru, J\f(au) ). 



The evolution equation (p4|) , 

^ h(t) = - n k (t) hit) , (73) 

appears then as the imaginary-time version of the Schrodinger equation with 
the self-adjoint, time-dependent, Hamiltonian Ti k . (The time-dependence of 
H,k comes from the ^-dependence of A/", k and a.) The formal "Dyson" 
solution of ([T3|) reads 

h(t) = Texp (- J^dt' W fc (0) MO)- (74) 

If the mean curvature k is positive, the Hamiltonian 7i k will be positive. 
If k is not positive, one can reduce the problem to a positive Hamiltonian 
TC' k = Hk + c (where c is a sufficiently large constant) by considering the 
modified variable h'(t) = e~ ct h{t). This positivity of the Hamiltonian makes 
it clear that the solution ( |74]) will exist for positive values of t under very 
mild conditions on the regularity of the initial data h(0 ) as a function of the 
surface coordinates £ . Like in the case of the heat equation, the evolution 
(|74|) for positive t will be well defined, and smoothing for all t . (The explicit 
solution of Eq. ( |T3"D given in the next subsection will show that h(t, £) is 
analytic in £.) 

In other words, we expect that our evolution, when considered in the 
outward direction (t > 0) will be well defined for all t , starting from a very 
general class of initial data (i.e. initial surfaces So ) with weak regularity 
properties. The surfaces S t generated by the outward evolution will be very 
smooth (probably analytic) even if the initial surface is barely regular. By 
contrast, the situation is completely different in the case of the evolution in 
the inward direction (t < 0). In this case, we expect that the initial data 
need in general to be an analytic function of £ , for the solution of ( [73]) to 
exist at all. The domain of existence of the solution (how "long" it exists for) 
will depend on how analytic the initial data is. This is precisely this type 
of feature that the family of norms entering the Nirenberg-Nishida theorem 



21, 22] are supposed to capture. 
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5.4 Explicit solution of the propagation equation for 
the perturbation of a monopole. 

We can actually write down directly (and explicitly) the solution h(t) of our 
propagation equation for the perturbation of a monopole, without even hav- 
ing to write this equation, in the case of the outward propagation (t > 0). 
To do so we consider, for some given value of t, the variation of the homoge- 
neous monopole Q t as described by the variation of both its inner (Eo) and 
outer (E t ) boundaries, expressed by the two height functions h(0) and h{t) , 
respectively. 

We recall that, in this variation, the total volume t spanned between 
the end surfaces is kept fixed. Therefore, the total mass M = \t\ of the 
homogeneous monopole (of unit density) does not vary, which implies that 
the variation of the external gravitational potential M/|x — x Q | must vanish. 
On the other hand, the variation of the external potential is equal to the 
potential generated by the variation of the volumic density of the monopole. 
The latter quantity is clearly made of two single layers: a layer of algebraic 
thickness + e h(t, £) on the outer boundary T, t and a layer of algebraic 
thickness — e h(0, £) on the internal boundary E . The problem is therefore, 
given the internal layer — eh(0,$,) , to determine the outer layer +eh(t, £) 
which is such that the external potential generated by these two layers vanish 
identically. This problem is solved by simple electrostatics considerations. 

Let us consider E t as a grounded conductor, and deposit on the surface 
E , enclosed within E t , a charge layer of surface density — e h(0, £) . By def- 
inition of a grounded conductor, the surface density induced on E t by the 
presence of the given layer distribution — e h(0, £) on E will be the unique 
solution +eh(t,£) of our problem. In mathematical terms, if Gs t (P,Q) de- 
notes the Dirichlet Green function of the domain interior to E t , the potential 
induced at any point P within E t by the layer — e h(0) on E reads 

v(P) = f G Et (P, Q) (-e h(0, Q)) rfS . (75) 
JQeS 

To simplify the notation we write h(0, Q) instead of a more explicit h(0, £q) . 

By definition of the Dirichlet Green function this potential vanishes when 
P G E t . The corresponding density induced on E t is given by the normal 
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derivative j^d n v on Et , so that, after division by e, 

h(t, P) = / rfE 4(P) /i(0, Q) . (76) 

The quantity cr^ (P) entering Eq. (|76|) is nothing but the previously defined 
surface density (|18D ( "harmonic measure" ) considered for E = E 4 and, as in 
subsection 3.3, with pole taken at an arbitrary point Q within the interior of 
E f (instead of the fixed origin O): 

4( p ) = -4^ <9nG Et (P, Q), (77) 
where <9 n denotes the outgoing normal derivative at P G E t . 



Eq. ( |76|) solves explicitly the problem of determining h{t) in terms of 
/i(0). The kernel K(t, P; 0, Q) defining the propagation of h(t), i.e. 
/i(t, P) = J K(t, P; 0, Q) /i(0, Q) c?Eo , written earlier in Eq. flT^| ) as a formal 
Dyson operator exponential, is in fact identical to cr^(P) . An alternative 
way of proving this result is to combine the "conservation" equation fl5T)|). 
which expresses the invariance of the gravitational energy of the monopole 
in an external potential U (x) , 



/ dZ t h{t,P) U{P) = [ rfS h(0,Q) U(Q) , (78) 
with the property (|2TD of the harmonic measure <t§(P) of determining the 



unique solution of the Dirichlet problem - i.e. the harmonic extension of 
data given on a surface E. More precisely, Eq. ([H]), written by changing 
O — > Q , allows us to express 

U (Q) = / rfE, a«(P) t/(P) , (79) 

if Q is in the interior of E 4 and U is harmonic within E t . Combining these 
two equations, and using the arbitrariness of the boundary data U (P) on E t , 
precisely yields the propagation equation ([76]) of h(t) . This argument based 
on the duality (^) between the perturbation h and the external potential 
U makes it clear why the kernel solving the (outward) propagation of h{t) 
(Eq. [76D is the adjoint of the kernel solving the (inward) harmonic extension 
problem (Eq. fT9|). 
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If the unperturbed family of surfaces St is analytic, the Green func- 
tion Gr£ t (P, Q) (which solves an elliptic problem) and its normal derivative, 
proportional to cr^ (P) , will be analytic in all their arguments. As a con- 
sequence, the outward propagated height h(t, P) , given by the integral ( [T6|) 
(where the points P and Q stay always separate when t ^ 0) will be, when 
t > 0, analytic in its arguments under very weak regularity assumptions for 
the initial height h(0, Q) 0. The kernel cr§ .(P) being positive, an initial dis- 
turbance which is everywhere positive on So is propagated to an everywhere 
positive h(t, P) on S t . Furthermore, even if the change h(0, Q) has compact 
support on S , the propagated change h(t, P) will, in general, be spread all 
over S 4 . Finally, the property ( |19|) of the kernel a ensures that the integral 
of h over S is preserved by the "time-evolution", which expresses that the 
volume of the monopole is kept constant (and also follows from taking U = 
const, in Eq. ([7Sp). 

The kernel representation ([H]) of the solution of the outward propagation 
problem does not help for solving the inward propagation one. By reversing 
the arguments above, it is clear that the inward propagation of a perturbation 
h(0) makes sense only if this perturbation is analytic in the spatial variables. 
In cases where, given an analytic h(0), the integral (|T6|) can be explicitly 
performed to give some well-defined analytic function h(t, for t > 0, its 
analytic continuation to negative values of t would, when it is possible, define 
the (unique) inward-propagated perturbation. 

Most of the results of this section hold, mutatis mutandis, when consid- 
ering the perturbation of a body which is gravitationally-equivalent to some 
fixed distribution pi (Q) . The perturbation may be caused, for example, by 
the variation of some inner boundary, and compensated by the appropriate 
variation of the outermost boundary. There is no difficulty in adapting the 
above reasoning to this CcLS6, clS well as to the perturbation of the more gen- 
eral (e.g. "Swiss-cheese-like" ) monopoles constructed at the end of subsection 

35 Note however that if h(0) has sharp variations in a given region of So (e.g. if it 
is angulous or even discontinuous) h(t) is expected to vary rapidly in the regions of St 
which are close to this region of So . Furthermore, h(t) may even present irregularities, 
reflecting those of h(0) , at points £ which are common to So and St . This situation is 
expected to occur, as discussed later in Section 6, whenever the surface So presents sharp 
angles (/3 < 7r/2) protruding towards its exterior. These points are expected to remain, at 
least for some time, fixed points of the successive surfaces S t built from S , both in the 
cases of outward and inward growths. 
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3.3. 



5.5 The perturbation of a spherical monopole. 

An illustrative example of the propagation of perturbation is the case of a 
spherically-symmetric monopole Q t , i.e. the volume contained between two 
concentric spheres of radii R(0) = Rq and R(t) . The center of gravity O 
must be the center of spherical symmetry. The parameter t is related to the 
radii by t — ^-(R 3 (t) — Rq) . If of the point P of S 4 at which we evaluate 
h(t), and Q those of the point Q at which the initial perturbation of S 
is located, the kernel solving the outward propagation problem (t > 0) is the 
Poisson kernel (^) with pole at Q, i.e. 

which tends asymptotically, for large t's, towards a uniform distribution 
a ~ 47r l R2( ^ , smoothing out all initial disturbances of S . 

An instructive way to see the effect of the propagation (especially in the 
inward direction for which it would be incorrect to use the kernel (|80|)) is 
to solve directly the "Schrodinger" equation (|73|) . For the case at hand, 
the harmonic measure is a^ t = (47r_R 2 ) -1 , the curvature k = 2/R and the 
Neumann operator J\f — t/R where i is the operator acting on functions on 
the sphere as does £ on the spherical harmonics Ye m (9,ip) : £Yg m = £Yg m . 
The explicit form of the evolution equation (1731) reads 



+ 2 K 



R 



R,9,<p) . (81) 



Its solution, given the "initial" perturbation of S described by the orthog- 
onal displacement h(Ro, 6, tp) , is as follows. If we decompose h in spherical 

harmonics, h(R,9,ip) = he m (R) Y^ m {6,ip) , we have 

e,m 

/R \ £+2 

W m {R) = h em (R ) . (82) 

Equation ([S2]) shows very explicitly the difference between the outward 
and the inward evolutions. In the outward case (R > Rq), the high-harmonic 
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fluctuations of the initial deformation are exponentially damped as 
exp [—(£ + 2) ln(R/Ro)} . This is a typical analyticizing transformation. In 
the inward case (R < Rq), the high-harmonics are enlarged by an exponen- 
tial factor exp [+(£ + 2) \n(R /R)} . Therefore the solution will exist (i.e. the 
spherical harmonics series will converge) only if the initial data were such 
that hi m (Ro) was an exponentially decreasing function of £, which means 
that the initial shape must be analytic. This illustrates the general properties 
of the solutions h(t) , that we have discussed in subsections 5.3 and 5.4. 

Evidently, one does not need the full apparatus of the Neumann operator, 
etc., to get Equation (RTF) and its solution Q8"2Tj. Let us briefly indicate how 



one can directly derive these results (generalized to any number of spatial 
dimensions N). Consider a system of polar coordinates (r,0) [where 6 = 
($x, . . . , 0jv-i) denote standard angular coordinates on the unit sphere in 
R N ] centered on the origin O which is to be the center of mass (and center of 
gravity). A would-be monopole having the topology of the domain contained 
between two nested spheres can be described by the polar equations of the 
inner (S ) and outer (Si) boundaries, namely r = po(0) and r = p\{0) , 
respectively. 

Let Yi m (6) denote a basis of spherical harmonics in N dimensions (with 
N — 2 "magnetic" quantum numbers m; in N = 2, m reduces to a ±1 
discrete index to distinguish e M from e~ %m , or cos £8 from sin £8). The 
multipole moments of the would-be monopole (with unit density) read 

Qim = 1^ d N * r l Y £m (0) ■ (83) 

J So 

Replacing d N jc — r N ~ l drdQg where dflg is the area element on the unit 
sphere in R N , and using the polar equations of the boundaries S and Si , 
we get 

Qim = J dtle Y em (6) jl^ [ Pl (0Y +N - p (Oy+ N ] . (84) 

In the case where So and Si are perturbed spheres of polar equations 
Pa(0) = Ra + th(R A , 6) (with A = 0, 1 , and £ < 1), one finds, to first 
order in e, and for £ > 1, 

Qim = e [ R^- 1 h^Ri) - R^ 1 h* em (R ) ] , (85) 
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where he m (RA) = J dftg Y^Jfi) ^(Ra, 0) is the £m spherical harmonics 
component of 1i(Ra, 0) . Equating Qg m (for £ > 1 ) to zero gives 



hi m (Ri) — (^j^- 



t+N-l 



(86) 



him (Roi) , 

which is the iV-dimensional generalization of ([52]). 

5.6 Towards a rigorous mathematical treatment of the 
growth of a monopole. 

To complete this Section dealing with the mathematical aspects of our con- 
struction, let us sketch the rewriting of the evolution equation ( p3|) describ- 
ing the growth of a monopole, as an evolution system that we anticipate 
to be amenable to a rigourous treatment along the assumptions of the ab- 
stract theorem of Nirenberg and Nishida. The basic idea is, like in many 
proofs of the Cauchy-Kovaleskaya theorem, to transform (E3) (whose right- 
hand side is nonlinear in the "dangerous" first-order terms, i.e. <9,f X and the 
pseudo differential dependence hidden in <te) into an equivalent "quasilinear" 
system. Essentially, such a system is obtained by introducing new indepen- 
dent variables, with corresponding evolution equations, until one obtains a 
closed, "quasilinear" system. We think that the system obtained by adding 
evolution equations for a, as well as for the ^-derivatives of the embedding 
map, p l a = d X l {t, £)/<9£ a , defines a "good quasilinear" system which should 



provide a basis for proving rigourously the existence of solutions of (|23|) in 
suitable functional spaces (of the type mentioned at the beginning of this 
Section). This system reads (suppressing the ^-dependence): 



dX l 
~dT 

da_ 
dt 



dpi 
dt 



a Tj{p> a ) 



A/* 



dp 

~d£« 



(* 2 ) , 



(87) 



da 
d^ 



n\p) + a 



dn l d p J /3 

dpj J¥ 
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The important point is that, here, n and k no longer denote the func- 
tional of the embedding (t, £ Q ) — > X l (t, £ Q ) defined above, but some explicit 
(rational) functions of and dp 3 [3 /d^ a obtained by replacing dX l /d^ a by 
p % a in the original expressions of the normal vector n [d X] and of the mean 
curvature k [d X, d 2 X] . Similarly, a no longer denotes the harmonic mea- 
sure, solution of an elliptic problem, but only some numerical function of £ . 
The non-locality (in £ ) of the system flB?D is all contained in the Neumann 
operator TVx associated with the surface x = X(t,£ a ). The information 
that a, p and dp/dt; are to be related to a° , dX/d£ and d 2 X/d^ 2 is 
fed in as suitable initial conditions taken at t = . The evolution system 
will then preserve this property 



The reasons why we are confident that one could establish the quasi- 
linearity of the system ([37]), i.e. the validity of inequalities of the type 
(|4"3"|), are the following. First, the facts that the normal n (p) does not 



depend on the derivatives of p with respect to and that the mean curva- 
ture k (p,dp/d^ a ) is linear in the derivatives dp/d^ a , ensure that all the 
algebraic terms on the right-hand sides of the system fl8T|) are quasi-linear. 
The only difficulty resides in the Neumann operator term. Here, the work 
of Verchota |^| should be sufficient to prove the quasilinearity of the term 



involving A/x ■ Indeed, if (using the notation of Ref. [p4}Q K denotes the 
in-surface double-layer potential operator and K* its adjoint, Verchota has 
shown that one could control the operators | + K and | + K* even when 
the domain Q t bounded by S t is Lipshitzian. Knowing that the inequali- 
ties he derived depend only of the Lipshitzian constant of Q t , and that M 

can be expressed in terms of Q + K^j or ( | + K*J and some domain- 
independent pseudo-differential operators, it seems clear to us that, with 
some work, one can establish the quasi- linearity of the term involving M . 



The system ( pTj) can also be applied to the evolution (|29[) of gravitationally- 
equivalent bodies. The information that a is to be related with the weighted 
harmonic measure (B8|) is to be fed in as an initial condition taken at t = 0. 



We wish to mention that the simple representation (^) of the linearized 
perturbation problem suggests an alternative route towards giving a rigourous 
proof of the existence of exact monopoles. If one rewrites the exact problem 
as the linearized problem (around the exact solution) plus nonlinear terms, it 
may be possible to use the form ( |T6T ) to set up a convergent iteration scheme 
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using at each step the exact "propagator" of the previous iteration. 

Finally, let us mention that Elliott and Janovsky have proven, for the 
case of two dimensions, the existence of a unique solution to an analogous 
problem - the injection of fluid in a Hele-Shaw cell - with, however, a differ- 
ent boundary condition near the origin [25|| . It would be interesting to see 
whether their method, based on an elliptic variational inequality reformula- 
tion of the moving boundary problem, can be adapted to our situation. 



5.7 The two-dimensional case. 

In the two-dimensional case one can prove the existence of particular classes 
of exact solutions of our nonlinear evolution problem by means of confor- 
mal mapping techniques. In N = 2 space dimensions, our problem is equiva- 
lent to the zero-surface-tension limit of the Saffman- Taylor problem pBl (see 
also the existence theorem of Ref. p5|). As shown, e.g. in Refs. 



some multi-parameter classes of exact solutions of the latter problem can 
be constructed by considering suitable t-dependent conformal mappings of 
the plane. Let us briefly discuss the application of this idea to our prob- 
lem (which has different boundary conditions than the usual Saffman- Taylor 
problem). 

One introduces the complex variables z = x + iy (where (x, y) are co- 
ordinates in the two-dimensional space) and C, = e~( G+lH \ where G(x,y) 
is the Green function, with pole at O, of the simply connected domain Q 
interior to some contour S, and where H(x,y) is its harmonic conjugate 
(dG/dx = dH/dy, dG/dy = — dH/dx). ( is a holomorphic function of 
z. The (Riemann) mapping z — > ( maps in a one-to-one conformal manner 
the domain Q onto the unit disc in the ( plane, the pole O being mapped 
on the center ( = of the disc (the logarithmic singularity of the Green 
function G ~ — In \ z — Zo\ in 2 dimensions ensures that ( ~ z — zo near 
the pole O). Actually, it is more convenient to consider the inverse mapping, 
z = f(Q , from the unit disc onto the considered domain Q. The t-evolution 
of a family of contours S t = dQ t can be described as the "dynamics" of 
t-dependent mappings z = ft(() , denoting the image of the unit disc 
under f t . 

The equation (E3l) governing the evolution of the boundary S 4 is equiv- 
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alent to enforcing the constraint 



Re[C9 c / t d t fn = lm[d e f t d t f t *} = ±- (88) 

I 71 

(independently of 9 and t) on the boundary of the unit disc, ( = e . The 
three ingredients needed to establish this constraint are: 

i) along the unit circle ( = e 10 , H is identical to — 9 ; 

ii) the harmonic measure is — tt- times the normal derivative of 
G, i.e. the tangential derivative of H; therefore the differential form ^- is 
mapped by ft into the harmonic measure a^ t dl, where dl is the line element 
along E t ; 

iii) the elementary area element swept dy dl when t is increased by dt 
is dA = dldh = dl a® t dt = dt ^ - also equal to ~ df t A dfl . This leads to 
Eq. (HD. 

Remarkably, one can find multi-parameter classes of mappings f(C',Pi, 
Pz,...,p n ) such that f t (() = f((;pi(t), . . . ,p n (t)) satisfies Eq. (H) when 
the Piit) satisfy some (complex) ordinary differential equations dpi/dt = 
Fi(pi, . . . ,p n ) p7|, ^|. The pj's parametrize singularities (zeroes and poles) 
of the derivative dft/dC,. In other words, one can reduce the dynamics of 
the contours S 4 to a particular dynamics of the critical points p^s. 

Let us illustrate this method by a very simple example. We consider a 
mapping of the form z = f t (() = a(t) C + | K*) C 2 • The derivative mapping 
df t /d( = a(t) + b(t) ( has a zero (and therefore a singularity) on the unit disc 
when \b(t)\ = \a(t)\. Let us consider the simple case where a(t) and b(t) are 
real. The constraint (|88|) is equivalent to the two differential equations 

1 d 1 
aa + - bo = — , 

2 2vr ' 

;j ( 89 ) 

ba + - ab = . 
2 

The second equation shows that a 2 b = C is a constant. The solution of the 
devolution is then implicitly given by 

1 o 1 , n lo 1 C t , , 

2 a +4 6 = 2 a = 2? +COmt ' (90) 
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The outward evolution of some initial "epicycloid" z = ao e %e + ~ 60 e 2ld (with 
I fro I < | fl o|) exhibits a growth of the "radius" a(t) and a fast decrease of the 
"ellipticity" b/a. The inward evolution develops a finger which becomes a 
cusp when \a\ = \b\ , i.e. in a finite time. This cusp is located away from the 
origin O, and locally described by y ~ (x — x ) 3 ^ 2 • We leave to future work a 
more thorough discussion of the usefulness of conformal mapping techniques 
in our context. 

6 Physical considerations. 

6.1 On the growth of aspherical monopoles. 

The examples of Section 4 and the results of Section 5 give a clear picture 
of the growth of a monopoleQ obtained by thickening, towards the outside, 
some initial surface S . It does not matter whether this one has a somewhat 
irregular shape (say a continuous shape with edges, corners, conical points, 
etc.). Generally speaking, the thickening of £ will, at first, be most impor- 
tant in the regions of £ located closest to the chosen origin O. The evolution 
stabilizes itself in that regions of fast growth recede fast away from O and 
thereby slow down their evolution. The shape of E becomes smoother and 
smoother and, in the long term, rounder and rounder. Note, however, that in 
cases where the initial surface E contains narrow gaps between sub-parts of 
E (e.g. when E is obtained by bending a "sausage- like" surface until the 
extremities get close to each other) - or if such a situation develops during 
the evolution - the subsequent outward evolution can develop a singularity 
( "collision" of surfaces) in a finite time. 

The fate of edges is interesting to consider. Consider on £ an edge 
with angle ft (opening towards the inside domain Q which includes O). 
The harmonic measure a is, initially, locally proportional to [|18[1 , 
where p denotes the distance (on £) away from the edge. This suggests that 
very acute edges, with ft < ~ , remain edgelike (with the same value of ft ) 
for a while, while the surface around curves out until, eventually, the edge 
disappears. Less acute edges, with ft > | (and even ft — |), are expected 
to disappear instantly (i.e. for any t > 0) because of the overlap (when 

36 Most of the following considerations also apply, mutatis mutandis, to the generalized 
case of the growth (03) of a "gravitationally-equivalent" body. 
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I < (3 < 7r) of the growing layers or (when (3 > it) because of their formally 
infinite growth rate. 

If we consider an infinite cylinder as initial X , with on the axis, we 
expect from Section 4 the development of an axially-symmetric bulge on the 
outside, with an exponentially-decreasing tail for large \z \ , proportionally to 

e~ 3 ~ L R' (the exact exponent ji in ( f4T|) being preserved during the evolution). 
During the thickening process the bulge should become rounder and rounder 
(and therefore less and less well approximated by the zeroth-order result of 
Section 4, represented earlier in Fig.l), tending asymptotically toward the 
shape of a sphere. 

On the other hand the growth of a monopole obtained by thickening, to- 
wards the inside, some initial surface £ corresponds to a completely different 
picture. Consider first the case where S is analytic (e.g. an ellipsoid or an 
infinite cylinder). Then the solution of Equation (|23|) should exist and stay 
analytic for some "time" t . However, there will be a tendency for bumps, 
and then "fingers", to grow towards the origin O. These fingers will become 
sharper and sharper. The process probably ends (assuming a "sufficiently 
analytic" initial shape £ ) by the development of an infinitely sharp spine at 
some distance of O, or, at the latest, when the fastest growing finger reaches 
the origin. In any case, there will be a finite "time" (i.e. volume) \t\ before 
the regular growth process breaks down, i.e. a maximum filling factor for 
generating a monopole with a given external surface. Even in simple cases 
(such as a sphere with O not being at the center, or an ellipsoid or a cylinder) 
it seems very difficult to estimate analytically this maximum filling factor. 

If the initial shape £ is non-analytic (and, say, contains corners, edges, 
conical points, . . .) it is not yet clear whether the inward growth process 
admits any rigourous solution. The formula given above for corners ( a oc 
p~P~ l ) suggests that very acute corners ( /3 < f ) will not fill, but will become 
more and more deeply incrusted within the solid body. On the other hand, a 
less acute corner ( (3 > |) might admit no uniquely defined inward evolution 
(when | < (3 < n one has the collision of two growing layers, while when 
(3 > 7r the tip is formally expected to grow at an infinite rate). 

In the case of small deviations from a spherical shape (or more gener- 
ally from a surface S), the results of subsections 5.4 and 5.5 show that the 
problem of inward propagation is closely akin to the problem of defining the 
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outward analytic extension of the inward analytic extension of a function 
given on the sphere (or on the surface £). As explained in subsection 5.4, 
the two problems are actually dual to each other. This duality leads to a 
simple explicit relation between the two problems in the case of deviations 
from a sphere. Indeed, if one starts from a function on the unit sphere, 

given, say, by its spherical harmonic expansion f(0) = J2 fe{9) > it admits 

j 

a unique inward harmonic extension within the sphere, f(r, 0) = ^ fl{&) 

e 

with r < 1. This well-defined inward harmonic extension is closely related^ 
to the result (^) (with R /R < 1 ) giving the well-defined outward growth 
of a nearly spherical initial £ . The problem of continuing the extension of 
f(r, 0) for r > 1 is thereby related to that of inwardly propagating a nearly 
spherical £ (Equation ( |8"2| ) with R /R> 1). 

Because of this analogy we might expect the inward growth problem to be 
mathematically ill-defined when the initial £ is non-analytic. Physically, the 
inward growth problem can make sense if one introduces new (regularizing) 
physical phenomena taking place at small distances. In fact, there is a large 
literature on such inward "Laplacian growth" problems (with cut off), and 
they lead to a complicated zoology of fractallike structure (for an entry into 
the literature see Ref. f58]] ). It does not seem that the resulting objects 
could be of practical interest for defining useful aspherical monopoles (or 
test masses). 

Let us finally remark that some interesting phenomena can take place in 
presence of symmetries (as for a parallepipedic or cubic box with the origin 
O at the center of symmetry). In such two or more bumps or fingers 

may grow at the same rate. In the case of inward growth, the symmetry 
between the fingers is unstable to small perturbations, which should lead to 
interesting phenomena of spontaneous symmetry breaking. 



37 More precisely the two problems can be identified if we consider the function 
f(R,9,<p) = {R/Ro) 2 h(R,9,<p) = E (R /R) e h £m {R ) Yi m (0,tp), and replace R /R 
by r. 
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6.2 On the practical realization of aspherical mono- 
poles. 

The construction we have presented could be used for actually fabricating 
aspherical (homogeneous) monopoles (or gravitationally similar test masses), 
in several different ways. 

A first type of approach^ consists in numerically computing beforehand 
an ideal shape, then machining the object according to this calculated shape, 
using some computer-driven machine. There are several practical methods 
one can use for computing an ideal shape. One is to follow strictly the 
differential construction given in Section 3, for instance by iterating an elec- 
trostatic code (solving the concerned Laplacian problem) in small successive 
"time" steps. Another method (a "spectral" method) is to parametrize the 
looked for final shape S t by a finite but large number of free parameters^], 
and to determine them by solving a correspondingly large number of con- 
straints, say the vanishing of as many multipole moments as possible. For 
instance, if we define the boundaries of the monopole by polar equations (as 
in subsection 5.5) the multipole moments are given by the simple formula 
fl81p (or more complicated ones if it is not enough to introduce two polar 
contours po(0,(f) and pi(9,<p)). If, for instance po(9,(p) is given, we can 
approximately solve for pi(8,(p) by expanding it as a linear combination of 
K basis functions (which may but do not need to be spherical harmonics) 
and numerically solve the first K — 1 equations Qe m = (£ > 1). The 
optimal choice of basis functions depends on the geometry of the problem 
and the expected qualitative shape of the solution. For instance the work of 
Section 4 on the cylinder suggests that one should use basis functions that 
capture the basic features of the lowest-order solution, namely a height func- 
tion above (or below) the surface p = R written as a function of ( = \z\/R 
which exhibits a smooth hump for |C| ~ 1 followed by an exponential decay 
for |C| > 1. 

38 This approach applies as well to the construction of gravitationally-equivalent bodies. 
By contrast the electrodeposition and brownian motion methods described below require 
more work to be applied to the case of Eq. (p9| ) . 

39 The "belted cylinder" approach introduced by the European STEP team |h| |, 
is a simple version of this approach. The better understanding we now have for the 
theoretically-optimal shapes can help to improve upon this approach (if nothing else, by 
considering suitable multi-belted cylinders with adequate tapering!). 
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A second type of approach for realizing a monopole may consist of using 
electrodeposition. Under certain conditions this may be a very direct way of 
realizing the inward thickening of an initial surface £ . Let us consider the 
space Q inside some given surface X to be (uniformly) filled with an elec- 
trolyte solution (of conductivity i] ) and let us apply a difference of potential 
between the cathode £ and an almost pointlike anode located at the origin 
O. The electric field E induces, in the quasi-stationary regime, a current 
density j = rj~E, satisfying V-j = except at the origin O. This current 
density obeys the same differential equation as E . Provided that we can 
consider the surface of the growing cathode as an equipotential, the magni- 
tude of j on E is proportional to the harmonic measure a . If we take for 
electrolyte a metallic solution, the thickness of the metal sheet deposited on 
the cathode £ per unit time, fixed by |j| , is proportional to a . Assuming 
the surface of the growing cathode to remain an equipotential under con- 
ditions allowing for the stable^ growth of a dense pattern, this method of 
electrodeposition can provide us with a practical mean to thicken £ precisely 
according to the law (|23|). 

An alternative way of getting this growth law is to have a stationary 
Brownian motion with a source of particles at O and to kill it when the 
particles first hit £ . The statistics of the hittings on £ is precisely given 



by the harmonic measure a (see e.g. Ref. [19|]), so that, if we conceive the 
particles as sticking on £ when they hit it, we have a Monte Carlo version 
of electrodeposition. This could be useful either as a computational method, 
or, possibly (if one can think of a suitable implementation), as a practical 
realization method. 



7 CONCLUSIONS. 

We have shown how to construct aspherical homogeneous bodies behaving 
exactly as point masses, both from the point of view of the gravitational field 

40 Many physical instabilities, in addition to the mathematical unstabilities mentioned 
above if the initial shape £ is not smooth enough, might alas show up. In particular, 
in order to avoid such instabilities we may have to ensure that the resistivity of the 
electrolyte is small enough, compared to that of the growing aggregate [29| , a requirement 
which tends to be in conflict with the other requirement that the surface of the growing 
aggregate remains an equipotential. 
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they generate, and from the point of view of their interaction with external 
gravitational fields. All their multipole moments (of order i > 1) vanish 
exactly. We have proven that these gravitational monopoles cannot occupy 
a "solid" region of space (in the sense of a solid sphere) but must possess (at 
least) one internal cavity. There is a very large flexibility in the construction 
of such objects (in any number of space dimensions). The free parameters 
entering this construction of monopoles in iV spatial dimensions (besides the 
value of the uniform density) are (i) an arbitrary "initial" hypersurface S 
bounding a connected open subset Q of R^; (ii) the arbitrary choice of the 
final center of mass (and "center of gravity" ) O within Q ; and (within some 
limits) the total volume of the constructed body. 

An extension of this method allows one to construct homogeneous bodies 
which are gravitationally equivalent (in the sense of having exactly the same 
infinite sequence of multipole moments) or gravitationally similar (with pro- 
portional sequences of multipoles) to any given body, or collection of bodies. 
There is a very large flexibility in the construction of these objects, given 
in particular the arbitrariness in the choice of the initial surface upon which 
the construction process is based. By combining and superposing the two 
constructions, we can obtain swiss-cheese-like monopoles having an arbitrary 
number of cavities, with arbitrarily chosen shapes. 

The basic construction consists of thickening the arbitrary initial hyper- 
surface X by depositing homogeneous layers with infinitesimal successive 
heights proportional to the harmonic measure, with respect to the fixed 
pole0 O, of the moving surface S t . The natural evolution parameter for 
this growing process is the (iV-dimensional) volume contained between X 
and S 4 . The topology of X is not restricted to be that of a sphere in R^; 
e.g. in 3 dimensions X can be a multi-handled "torus". The growth pro- 
cess has very different properties when performed outwardly, or inwardlyf 2 ], 
with respect to X . Our study of the linearized evolution equation (a kind 
of imaginary-time Schrodinger equation, i.e. similar to the heat equation), 
and of its solutions, shows that the outward growth process is stable and 

41 The construction of gravitationally-equivalent bodies is similar, provided one replaces 
the harmonic measure cr£ t by the weighted measure u^t of Eq. (j^) . 

42 Though we phrased our construction in the text as being either outward or inward, it 
is evident that one may want, in some practical applications, to thicken some "initial" X 
on both sides at once. 
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smoothing. It is probably well defined even if the initial hypersurface has 
a very limited regularity. By contrast, the inward growth process is, proba- 
bly, mathematically well defined only for analytic initial data. We presented 
the program of a rigourous mathematical proof of the existence of aspherical 
monopoles without, however, completing such a proof. We also exhibited 
explicitly a simple example of a two-dimensional aspherical monopole. 

We indicated several practical ways of realizing gravitational monopoles 
(including computer-aided machining and electrodeposition). In particular, 
once one knows, from the general construction method, that such objects ex- 
ist with many possible shapes, one can try to find numerical approximations 
to some of them (as Barrett did some years ago) without necessarily fol- 
lowing into the footsteps of the harmonic measure construction (see Section 
6.2). Though our method generally assumes that the domain Q is bounded, 
it can also generate unbounded monopoles, with cylindriclike internal cavi- 
ties, and exponentially-decreasing thicknesses at infinity. 

Apart from its conceptual interest, our construction might be of practical 
interest for optimizing the shape of test masses in high-precision Equivalence 
Principle experiments, such as STEP. By suppressing the coupling to gravity 
gradients, it allows one to design, with great flexibility in the choice of shapes, 
differential accelerometers, made of nested bodies, which are insensitive to 
external gravitational disturbances. Even very close disturbing masses would 
have no effect^ on the differential accelerometer, which is then "gravitation- 
ally screened" from all external gravity gradients. Truncated versions of the 
cylindriclike monopoles with exponentially falling off thicknesses may define 
particularly useful shapes for the STEP experiment. If their length L is large 
enough compared to their radius R such cylinders would be exponentially 
insensitive to disturbing masses^, even if these are very close. 

Using the extended construction method described in subsection 3.3 we 
can choose any (non-monopole) shape for some inner body, e.g a straight 
cylinder, and construct an enclosing (or cylinderlike) outer body which has 
exactly (or almost exactly) the same set of reduced multipole moments 

43 When the monopoles are concentric. 

44 Our lowest-order solution - and its successive iterations - indicate that the quality 
factor of the decoupling goes like e 2A05L / 2R ( increasing by about an order of magnitude 
each time the total length L is increased by 2R . 
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Qe /M as the inner body. This is enough to ensure that the differential 
accelerometer they define is insensitive to all external gravity gradients. 

Finally our techniques may also allow one to devise multi-nested differen- 
tial accelerometers, made of "Russian doll" type masses, (almost) insensitive 
to external gravitational disturbances. This could be important for conceiv- 
ing and implementing more compact (and cheaper) versions of the STEP 
concept. Altogether, the optimization of the test masses used in differential 
accelerometers should allow for a significant increase of their performances, 
as required for very precise tests of the Equivalence Principle. 
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Figure caption. 

Fig. 1: Surface density for a thin cylindrical gravitational monopole of radius 
R, as a function of the longitudinal coordinate z. The function F (z/R) , de- 
fined by Eqs. ( j37"D or (f40|) , determines the mass density a (z) = F (z/R) 
All (/ > 1) multipole moments vanish identically. 86 % of the mass corre- 
spond to \z\< R, 98.7 % to \z \ <2R. 
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